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Foreword 


Knowledge  of  the  space  and  time  evolution  of  environmental  parameters 
in  the  ocean  can  significantly  improve  acoustic  surveillance  capabilities, 
performance  of  weapons  systems,  search  and  rescue  planning,  and  other  aspects 
of  naval  operations.  Observations  alone  cannot  provide  this  information  and 
must  be  supplemented  by  the  use  of  numerical  forecast  models.  The  Naval 
Ocean  Research  and  Development  Activity  is  designing  and  developing  a  hier¬ 
archy  of  numerical  forecast  models,  including  multilayer  hydrodynamic  models, 
for  Navy  use  in  predicting  the  ocean  environment.  The  energy  analysis  of 
forecast  results  is  a  major  diagnostic  technique  in  testing  forecast  results  and 
must  be  formulated  specifically  for  each  numerical  forecast  model. 


A.  C.  Esau,  Captain,  USN 
Commanding  Officer,  NORDA 


Executive  summary 


Energy  analysis  is  a  common  technique  for  examining  the  physical  behavior 
of  fluid  systems.  The  specifications  of  energies  and  their  time  variations  are 
not  unique  but  depend  on  the  dynamical  model  being  studied.  A  dynamical 
model  consists  of  a  set  of  equations  giving  the  mathematical  representation 
of  the  relevant  physical  laws,  together  with  the  particular  assumptions  and 
approximations  being  applied.  Here  two  dynamical  models,  useful  in  studying 
ocean  circulations,  are  treated:  the  inviscid  multilayer  hydrodynamic  model 
and  the  inviscid  multilayer  reduced-gravy  hydrodynamic  model  or,  for  brevity, 
the  finite  depth  model  and  the  reduced-gravity  model.  The  two  models  are 
identical  except  that  the  velocity  is  required  to  be  zero  in  the  bottom  layer 
of  the  reduced-gravity  model.  The  reduced-gravity  model  retains  most  of  the 
physical  properties  of  the  finite  depth  model  but  has  considerably  fewer 
computational  requirements. 

Two  different  formulations  of  the  energetics  are  obtained:  layer-layer  and 
layer-interface.  In  both,  the  kinetic  energy  per  unit  area  of  a  layer  is  specified. 
In  the  layer-layer  formulation  the  available  potential  energy  per  unit  area  for 
a  layer  is  specified  while  the  layer-interface  formulation  uses  the  available 
potential  energy  per  unit  area  associated  with  an  interface.  In  addition,  the 
types  of  energy  transfer  possible  are  different  for  the  two  formulations.  For 
the  layer-layer  formulation  kinetic  energy  and  available  potential  energy  are 
exchanged  within  a  layer  and  available  potential  energy  is  exchanged  between 
adjacent  layers.  There  is  no  transfer  of  kinetic  energy  between  layers.  For 
the  layer-interface  formulation  energy  is  transferred  between  the  kinetic  energy 
in  a  layer  and  the  available  potential  energy  in  an  adjacent  interface,  above 
or  below  the  layer.  In  each  formulation  the  magnitudes  of  the  energy  transfer 
terms  differ  for  the  two  models.  In  addition,  for  the  layer-interface  formulation 
the  reduced-gravity  model,  but  not  the  finite  depth  model,  permits  the  exchange 
of  available  potential  energy  between  consecutive  interfaces.  This  latter 
exchange  is  a  consequence  of  the  required  coupling  of  interface  deviations 
in  reduced-gravity  models.  Finally,  both  formulations  are  valid  mathematically 
and  physically  so  the  selection  of  one  depends  on  what  types  of  energy  transfers 
the  user  wishes  to  examine. 

The  energetics,  or  energy  analysis,  is  presented  for  a  general  multilayer  case 
for  both  models  and  both  formulations.  Also,  the  energetics  in  both 
formulations  are  considered  for  specific  cases:  the  one-,  two-,  and  three-layer 
cases  for  the  finite  depth  model  and  the  one  and  two  active  layer  cases  for 
the  reduced-gravity  model. 
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Energy  Transfers  for  Multilayer  Hydrodynamic  Ocean  Models: 
General  Formulation  and  Specific  Examples 


Glossary 


i  available  potential  energy  per  unit  area  in  layer  i 
Aj  available  potential  energy  associated  with  the  interface  above  layer  i 

d/dt  3/9t  +  d/dx  +  d/dy,  total  time  derivative 

energy  flux  divergence  for  layer  i 
f  Coriolis  parameter 

g  acceleration  due  to  gravity 

thickness  of  layer  i  when  system  is  at  rest 
hj  actual  layer  thickens 

interface  deviation  at  top  of  layer  i 
i,j  subscripts  referring  to  a  particular  layer  or  interface 
i,j  unit  vectors  in  the  x,  y  directions 

kinetic  energy  per  unit  area  for  layer  i 
unit  vector  normal  to  side  walls  in  layer  i 


N  total  number  of  layers  or  interfaces 
depth-dependent  pressure  for  layer  i 
Pj  depth-independent  pressure  for  layer  i 

p^  modified  depth-independent  pressure  for  layer  i,  used  in  reduced-gravity  model 
potential  energy  per  unit  area  for  layer  i 
minimum  potential  energy  per  unit  area  for  layer  i 

Pj  density  in  layer  i 

RHS  right  hand  side 


1 


S  area  of  side  walls  in  layer  i 

1 

t  time 

V  =  ui  +  vj 

w .  r .  t .  with  respect  to 
X  east-west  distance,  positive  to  the  east 

y  north-south  distance,  positive  to  the  north 

z  vertical  distance,  positive  upward  and  zero  at  the  area-averaged  ocean  top 

z.  actual  depth  at  bottom  of  layer  i 

1 

^  =  i  a/ax  -I-  j  a/ay 

3/3t  local  time  derivative 

^  I  angle  brackets  indicate  an  average  over  the  area  of  the  system 
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PART  A;  INTRODUCTION 


Numerical  forecasting  of  the  ocean  environment  requires  first  a 
consistent  set  of  equations  that  express  the  relevant  physical  laws  for  the 
dynamics  of  geophysical  fluids  together  with  appropriate  initial  and  boundary 
conditions  and  then  a  rapid  and  accurate  numerical  method  for  integrating  the 
equations  in  space  and  time.  The  physical  laws  for  the  ocean,  in  general, 
include  conservation  laws  for  mass,  momentum,  and  salinity,  an  equation  of 
state  and  a  thermodynamic  relation.  A  unique  mathematical  representation  of 
these  physical  laws  does  not  exist.  Furthermore,  due  to  the  complexity  of 
these  equations,  various  assumptions  and  approximations  may  be  introduced 
which  simplify  the  equations  or  even  reduce  the  number  of  equations.  The 
numerical  integration  method  is  then  tailored  to  the  chosen  set  of  equations. 
This  set  of  equations  together  with  the  particular  integration  method  is 
referred  to  as  a  numerical  forecast  model.  Any  such  model  may  be  subject  to 
different  sets  of  boundary  and  initial  conditions  provided  each  set  used  is 
consistent  with  the  given  model. 

The  validity  of  the  numerical  forecast  generated  by  a  particular 
numerical  model  with  a  selected  set  of  boundary  and  initial  conditions  is 
judged  in  part  by  comparison  with  observed  conditions  in  the  ocean  region. 
Improving  the  forecast  requires  changing  the  boundary  and  initial  conditions 
applied  or  changing  the  numerical  model  itself.  Such  changes  must  be  based  on 
a  thorough  understanding  of  the  dominant  physical  mechanisms  incorporated  in 
the  model  and  in  the  applied  boundary  and  initial  conditions.  Analysis  of  the 
energetics  in  the  forecast  region  is  one  technique  used  to  fathom  the  physical 
processes  present. 

An  introduction  to  energy  analysis  in  geophysical  fluid  dynamics  can  be 
found  in  a  number  of  textbooks  (e.g.,  Dutton,  1976;  Gill,  1982;  Holton,  1972; 
Pedlosky,1987;  and  Phillips,  1977).  The  three  types  of  energy  considered  are 
kinetic,  potential,  and  internal,  but  the  specific  formulation  of  each  is 
dependent  on  the  given  set  of  equations  being  used.  Thus,  the  energetics 
analysis  for  a  numerical  model  is  dependent  on  the  chosen  set  of  equations  and 
must  be  rederived  for  any  change  in  that  set.  One  of  the  objectives  of  the 
Naval  Ocean  Research  and  Development  Activity's  (NORDA)  FY  86-87  Tactical 
Oceanography  Project  is  to  develop  a  general  energetics  analysis  for  the 
multilayer  hydrodynamic  numerical  model  developed  originally  by  Hurlburt  and 
Thompson  (1980),  hereafter  referred  to  as  H  &  T. 

The  development  of  a  general  energy  analysis  for  a  given  model  begins 
with  the  formulation  of  the  types  of  energy  present.  Next,  the  change  of 
these  energies  with  respect  to  time  (w.r.t.)  is  obtained.  This  step  usually 
requires  considerable  mathematical  manipulation  to  obtain  a  physically 
reasonable  formulation  for  the  transfer  of  energy  between  the  various  types  of 
energy  present  and  between  the  various  regions  in  the  fluid  domain.  Finally, 
the  energy  analysis  for  numerical  models  usually  involves  some  type  of  space 
and/or  time  averaging.  This  report  includes  only  the  first  two  steps  which 
are  interrelated  and  lengthy;  the  averaging  techniques  will  be  presented  in  a 
later  report. 

The  set  of  equations  appropriate  for  an  inviscid  multilayer  hydrodynamic 
model  is  derived  in  the  Appendix.  These  equations,  valid  for  all  layers,  form 
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the  basis  for  the  general  energy  analysis  and  are  listed  in  Table  1  for  ready 
reference.  For  some  research  investigations,  an  additional  assumption  is 
made:  namely,  that  the  bottom  layer,  N,  is  assumed  to  be  of  infinite  depth. 
The  numerical  model  with  this  additional  assumption  is  called  an  inviscid 
multilayer  reduced-gravity  hydrodynamic  model.  For  brevity,  the  original 
model  is  called  the  finite  depth  model  and  the  second  model  is  called  the 
reduced-gravity  model.  The  energy  analysis  for  both  models  is  developed,  but 
a  derivation  of  the  changes  in  the  basic  set  of  equations  for  the  reduced- 
gravity  model  is  delayed  until  section  B.5. 


Table  1.  Basic  equations  for  the  finite  depth  model 


3u. 
_ 1 

at 


+ 


fv. 


1 


(1) 

(2) 

(3) 

(4) 


The  set  of  equations  in  Table  1  differ  from  those  found  in  H  &  T  in 
several  respects.  First,  the  equations  in  H  &  T  are  in  transport  form  and 
those  in  Table  1  are  in  momentum  form.  The  momentum  form,  rather  than  the 
transport  form,  is  used  in  this  report  because  velocity,  not  transport,  is 
used  in  formulating  kinetic  energy.  Second,  the  equations  in  Table  1  do  not 
contain  any  external  forcing,  such  as  winds,  nor  do  they  have  any  viscous 
effects  or  dissipation  included.  Such  sources  and  sinks  determine  the 
resulting  circulation  in  hydrodynamic  models  but  their  inclusion  in  any  energy 
analysis  is  straightforward:  the  sources  increase  the  kinetic  energy  and  the 
sinks  decrease  the  kinetic  energy.  The  difficult  part  of  any  energy  analysis 
is  determining  a  mathematically  consistent  and  physically  reasonable 
formulation  for  the  tranfers  of  energy.  Once  this  formulation  is  determined, 
adding  the  source  and  sink  terms  is  quite  simple.  Therefore,  the  source  and 
sink  terms  are  omitted  in  this  report  to  simplify  the  mathematical  analysis 
included. 

The  energy  analysis  for  layered  models  is  developed  in  section  B.l 
without  reference  to  the  number  of  layers  present.  Various  ways  of  treating 
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The  energy  analysis  for  layered  models  is  developed  in  section  B.l 
without  reference  to  the  number  of  layers  present.  Various  ways  of  treating 
the  kinetic  energy  and  the  available  potential  energy  are  presented  in  the 
following  sections  in  Part  B  and  equations  showing  the  change  of  these 
energies  w.r.t.  time  are  derived.  One  significant  result  of  this 
investigation  is  that  there  are  a  variety  of  ways  to  express  the  change  w.r.t. 
time  of  both  the  kinetic  and  available  potential  energies.  Two  formulations 
are  selected  for  further  consideration:  the  first  or  layer-layer  uses  the 
available  potential  energy  in  a  layer  while  the  second  or  layer-interface 
treats  available  potential  energy  associated  with  each  interface.  In  both 
treatments  the  kinetic  energy  for  a  layer  is  used. 

Many  forecasts  made  with  multilayer  models  actually  use  only  one  or  two 
layers.  Also  the  energy  relations  for  the  top  (bottom)  layer  differ  from 
those  for  the  internal  layers  because  there  is  no  layer  above  (below)  with 
which  to  exchange  energy.  Therefore,  in  Parts  C  (layer-layer  formulation) 
and  D  (layer-interface  formulation)  the  specific  energy  relations  for  one, ^ 
two,  and  three  layers  are  given  for  the  finite  depth  model  and  for  one  active 
and  two  active  layers  for  the  reduced-gravity  model.  Each  section  in  Parts 
C  and  D  is  designed  to  stand  alone;  therefore,  there  is  some  overlap  in  the 
descriptions  presented  in  each  section. 

This  report  is  designed  for  two  purposes:  first,  to  delineate 
the  energy  analysis  formulations  best  suited  to  the  inviscid  multilayer 
hydrodynamic  model  and  to  the  inviscid  multilayer  reduced-gravity  hydrodynamic 
model;  second,  to  serve  as  a  training  manual  for  those  scientists  who  are 
unfamiliar  with  the  techniques  of  energy  analysis  used  with  forecast  models. 
The  second  purpose  dictates  that  a  considerable  amount  of  explanation  and 
mathematical  detail  be  included  that  generally  is  eliminated  for  brevity. 
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Part  B,  Section  1 


Kinetic,  potential,  and  internal  are  the  three  types  of  energy  commonly 
considered  in  the  dynamics  of  geophysical  fluids.  In  hydrodynamic  models,  the 
temperature  is  not  a  variable;  therefore,  internal  energy  is  not  considered. 
Furthermore,  the  models  being  analyzed  are  layer  models  so  it  seems  reasonable 
to  examine  kinetic  energy  and  potential  energy  in  a  layer,  rather  than  for  the 
entire  volume  of  the  ocean  region.  The  structure  of  the  layers  in  the  finite 
depth  model  is  shown  in  Figure  1. 

The  kinetic  energy  density  in  layer  i  is  defined  as  p.  fu.^  +  v.^1/2. 

The  column  integrated  kinetic  energy  density  or  the  kinetic  energy  per  unit 
area  in  a  layer  i  is  K^  where 


K. 

1 


^  2  ‘^iPi 


+  V 


(1.1) 


All  symbols  are  defined  in  the  Glossary. 


The 

derived 


change  w.r.t.  time 
in  the  Appendix  and 


of  K.  can  be  derived  from  the  basic  equations 
listed  in  Table  1; 


fau  ) 

pUi^ 

rau.) 

1 

'ap.^ 

U  J 

+  u. 

1 

.3x  J 

+  V. 

1 

J. 

Uy  J 

-  fv.  =  -  — 

1  p. 

lax  J 

rav.) 

1 

rav.) 

rav.) 

1 

1 

lat  J 

+  U. 

1 

lax  J 

+  V  . 

1 

lay  J 

rah.) 

rah.) 

rah.) 

fau.) 

fav.) 

1 

1 

1 

1 

1 

.M 

CO 

+  U. 

1 

lax  J 

'^ilay  J 

+  h . 

1 

lax  J 

+  h , 

1 

lay  J 

i 

p.  .  p.g(Pj  -  z)  -  (Pj  -  p.)Bllj 


(1.2) 


(1.3) 


(1.4) 

(1.5) 


Multiply  (1.2)  by  h.p.u.  and  (1.3)  by  h.p.v..  Add  the  resulting  equations  to 
obtain  ^  ^  ^  ^  ^  ^ 


N  It  [i  *  ’i^)]  *  "i 


1(2  2) 

Next  multiply  (1.4)  by  j  p^  lu^  +  v^  I  to  yield 


(1.6) 


(1.7) 
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Add  (1.6)  to  (1.7)  to  obtain 


aK. 

1 

at 


(1.8) 


or 


aK. 

1 

at 


^•v.K.  -  ^‘V.h.p.  +  p.^'v.h.  . 
11  1  1*^1  ^1  11 

(1)  (2)  (3) 


(1.9) 


Term  1  is  the  layer-integrated  flux  of  kinetic  energy  or  the  layer-integrated 
kinetic  energy  flux  density.  Term  2  is  the  layer-integrated  pressure  work 
term.  Sometimes  terms  1  and  2  are  combined  to  form  the  layer-integrated 
energy  flux  divergence.  If  the  layer-integrated  energy  flux  divergence  is 
integrated  over  the  horizontal  area,  it  represents  the  flux  of  energy  through 
the  side  walls.  Term  3  is  discussed  later. 

It  may  appear  that  (1.9)  contains  a  dependence  on  the  depth,  z,  which 
(1.8)  does  not.  It  is  readily  shown  that  (1.9)  does  not  depend  on  z. 

From  (A. 7) 


P^(x,y,z,t)  =  p.(x,y,t)  -  p.gz  . 


Substituting  (1.10)  into  (1.9)  gives 


SK, 

3t 


i  =  -  ^-v-K.  -  ^’V-h-p.  +  p.^'V.h.  , 


1  1 


1  i‘  1  ‘1  11 


(1.10) 


(1.11) 


which  has  no  z  dependence.  Therefore,  it  does  not  matter  whether  p.  or  p.  is 
used  in  the  energy  calculations;  there  is  no  z  dependence  in  K..  However,  p. 
will  be  used  in  the  remainder  of  the  report  to  simplify  the  calculations.  ^ 

The  potential  energy  density  in  layer  i  is  defined  as  p^gz.  The 
potential  energy  per  unit  area  in  a  layer  is  defined  as 


P. 

1 


p^gz  dz 


(1.12) 


or 
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(1.13) 


P. 

1 


where  is  the  depth  at  the  bottom  of  layer  i.  Note  that  z^<0  for  i>l. 
From  Figure  1,  can  be  written  as 


2  . 
1 


=  h 


i  +  l 


1 

H. 

J=1  3 


for  i>l 


(1.14) 


and 


z 

o 


(1.15) 


For  fluids  most  of  the  potential  energy  is  not  available  for  conversion  to 
other  forms  of  energy.  The  system  shown  in  Figure  1  is  stably  stratified;  the 
density  in  a  given  layer  is  less  than  in  the  layer  below.  Therefore,  the 

minimum  potential  energy,  P.,  occurs  when  the  X],'s  are  all  zero  or,  from 
(1.14)  and  (1.15)  ^  ^ 


'■i '  i  "is  ^ 

where 


^z.  ^  =  -  L  H.  for  i>l  (1.17) 

and 


^z  1=0. 

^  o  ^ 


(1.18) 


Note  that  both  P.  and  P.  are  less  than  zero.  This  result  occurs  because  zsO 
at  the  top,  not  at  the  feottom  of  the  fluid.  It  is  not  convenient  to  set  z=0 
at  the  bottom  of  the  fluid  because  the  bottom  layer  is  infinitely  deep  for  the 
reduced-gravity  model  to  be  considered  in  section  B.5. 

The  amount  of  potential  energy  that  is  available  to  be  converted  to 
kinetic  energy  is  the  difference  between  the  actual  potential  energy  and  the 
minimum  potential  energy.  Therefore,  the  available  potential  energy  per  unit 
area  in  a  layer  is 
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(1.19) 


A.  E  P.  -  P. 
Ill 


Substituting  (1.13)  and  (1.16)  into  (1.19)  gives 

\  * 
Using  (1.14),  (1.15),  (1.17)  and  (1.18)  in  (1.20)  yields 

(- 


A. 

1 


fPig 


r 

Y  i-1  ] 

2 

i-1  i 

2^ 

n.  -  EH. 
^  i-1  ^ 

A  J  -  J-  J 

_ 

-EH. 

j=l 

. 

2  2  ) 
\  (  \  '  T 


n.  n  -  I  H. 
j=l 


-EH. 

j=l  ' 


i>l 


(1.20) 


(1.21) 


or 


A.  = 
1 


fPiSl 


2 


1 

n.  -  2n.  2  H. 

V  "  "  j=i  ^ 


'i+l 


1 

E 


H. 


2ri.  1  -  - 

'"1  j.i  j 


For  the  top  layer,  i=l,  (1.22)  becomes 

*1  =  (¥]  K'  -  h'  -  ^”2  "i)] 


i=l  . 


i>l  .  (1.22) 


(1.23) 


The  change  of  available  potential  energy  w.r.t.  time  is  found  by  taking 
9/3t  of  (1.22)  and  (1.23)  to  get 


3A. 

W  "  '’i^ 


i-1  )  30.  ^ 

1 


0.  -  r  H. 

j=i  ' 


3t 


i+l 


1 

-EH. 

i=i  ' 


an 


i+l 


3t 


i>l  (1.24) 


and 

3A^ 

3t“  "  ^1® 


30i  .  v902 

.'^l  JT  ~  K  "  ^iJsF 


i=l  . 


(1.25) 
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The  total  energy  per  unit  area  in  a  layer  is  defined  as  the  sum  of  the 
kinetic  energy  per  unit  area  and  the  available  potential  energy  per  unit  area 
or  K^+Aj^.  Adding  (1.11)  and  (1.24)  gives  the  change  w.r.t.  time  of  the  total 
energy  per  unit  area  in  a  layer  as 


■It  [k.  +  A.l  =  -  ^'V.  fk.  +  h.p.l  +  p.^-v.h. 
9t  (  1  i)  ill  i^iJ  *^1  11 

an, 


+  pjg 


f  i-1 

n.  -  r  H. 
h  j=i  ^ 


at 


^.1 


1 

z 


j=i 


H. 

3 


an.  1 ' 

1+1 

at 


(1.26) 


The  first  term  on  the  RHS  of  (1.26)  represents  the  horizontal  transport  of 
energy  into  the  area.  There  are  no  external  sources  and  sinks  of  energy  for 
the  finite  depth  model  described  in  the  Appendix.  Therefore,  the  remaining 
terms  on  the  RHS  of  (1.26)  must  represent  vertical  transport  of  energy  or  the 
exchange  of  energy  between  layers  in  the  system.  However,  it  is  not 
immediately  obvious  from  (1.26)  which  terms  represent  transfers  between 
layers.  Furthermore,  transfers  between  kinetic  and  available  potential 
energies  within  a  layer  should  be  possible,  but  these  transfers  are  not 
obvious  from  (1.11)  and  (1.24).  The  remaining  sections  in  Part  B  examine 
various  ways  to  formulate  energy  transfers  between  layers  and  within  a  layer 
for  both  the  finite  depth  and  the  reduced-gravity  models. 
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Part  B,  Section  2 


In  this  section  the  specific  equations  for  energy  transfer  between  layers 
and  within  a  layer  are  obtained  in  a  layer-layer  formulation  for  a  multilayer 
model. 


From  section  B.l,  (1.11)  and  (1.24),  the  equations  for  the  local  change 
of  kinetic  energy  and  of  available  potential  energy  with  w.r.t.  are 


aK. 

1 

at 


+  h .  p . 


(2.1) 


and 


aA. 

1 

at 


Pig 


^i 


i-l 

E 

j=l 


H. 

J 


Figure  1  shows  that  the  relation  of  h^  to  1)^ 


and  H.  is 
1 


(2.2) 


h.  =  n.  +  H.  -  t).  1  l<i<N  ,  (2.3) 

111  1+1  -  ’  ^  ' 

while 


(2.4) 


Substituting  (1.10)  into  (1.5)  gives 


Pi  .  p.eiii  -  (p.  -p.)gh.  , 

while  (1.4)  is 


(2.5) 


(2.6) 


The  preceding  set  of  equations,  (2.1)  -  (2.6),  can  be  reduced  by  eliminating 
the  h's.  Sura  (2.3)  over  all  N-1  layers  and  add  (2.4)  to  obtain 

N  N 

E  h.  =  hi  +  2  H.  (2.7) 

j=l  ^  ^  j=l  ^ 


11 


or 


N 


1 

I 

j=l 


,  -  H.  . 

(2.8) 

)  3J 

summed  from  layer  1 

to  layer  i  where  i<N,  the  results  are 

i 

-  n.  .  +  z  H.  . 
j=i  J 

(2.9) 

by  using  (2.8)  in 

(2.9)  to  find 

r  -s  ^  f 

h.  -  H.  +  Z  H. 

-  h.  =  Z  h.  -  H  . 

(2.10) 

V  3  3)  j=i  J 

j-i.i 

’i  -  h  -  “3]  • 


Use  (2.10)  in  (2.3)  to  obtain 

N 
I 

j  =  i 

Substituting  (2.8)  in  (2.5)  gives 

N  N  i-1 

p.  =  p.g  Z  h.-  p.g  L  H.  +  Z  P'gh.  . 

j=i  ^  j=l  ^  j=l  ^  ^ 

Use  (2.11)  and  (2.10)  in  (2.2)  to  get 


8A. 

W  "  ^i® 


Z  h 

.  -  H.l 

LI  j=i^' 

]  3) 

i-1 
Z  H. 

j=l  ^ 


^  .  r  N  'i 

d 

Z  h . 

9t 

/ 

.  j=i 

f  ^  ^ 

i  ) 

a 

'  N  y 

Z  h.  -  H. 

-  EH. 

0 

Z  h. 

i=i  'J 

Jt 

,  j=i+l 

or 


(2.11) 


(2.12) 


(2.13) 
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9A. 
_ 1 

3t 


'  N  N  ' 

N  ) 

=  Pig 

Eh.-  EH. 

i-i  J  i-1  ^ 

A  3=1  3-1  > 

0 

It 

E  h. 

.  j=i 

'  N 

N  ) 

3 

'  N  y 

E  h. 

-EH. 

0 

E  h. 

w  j=i+l  ^ 

3=1 

9t 

,  j=i+l 

(2.14) 


The  terms  in  2.14  may  be  rearranged  to  yield 


3A. 

"St”  ^ 


N  N 

E  h.  -  EH. 
J 


j=i+l 


i=l 


3h .  N  3h  .  ■ 

sT  ^  ^i.f .  3t^ 


(2.15) 


All  terms  on  the  RHS  of  (2.15)  involve  time  derivatives  of  h^.  Multiply  (2.6) 
by  and  substitute  in  (2.1)  to  find 


3K.  ^  ^  ^ 

— -i  =  -  V'V.  K.  +  h.p.  -  p.  -5-;—  . 

3t  1(1  11)  1  3t 


(2.16) 


3h. 

1 


Add  and  subtract  the  term  p^  qY"  from  (2.15)  to  find 


3A. 
_ 1 

3t 


= 


3h. 

3h. 

f  ^ 

N  '' 

3h, 

N 

3h.' 

1 

3t 

Pig 

E  h.  ■ 
. .  j=i+l  ^ 

-  E  H. 
3=1 

1 

IF 

+  h.  E 

1 .  . 

3  =  1 

3 

3t 

(2.17) 


Replace  the  second  p^on  the  RHS  of  (2.17)  with  (2.12)  to  yield 


3A. 

1 


3h. 


3t  ^i  3t 


'  N  N  ^ 

i-1 

Eh.-  EH. 

+  ^  PnSb. 

V  j=i  ^  3=1 

.  .  1  1 
3=1  -^J 

3h. 
_ 1 

3t 


+  P^g 


■f  N 

N  ^ 

E  h.  - 

E  H. 

A  j=i+l  ^ 

j=l 

3h, 


N  3h. 


3=1 


(2.18) 


or 


3A. 


3h. 


N  3h. 


3t 


-  '  fi  ar  *  ^  - 

J=i+1 


i-1 

E  p.gh. 
•  •  3  3 

3  =  1 


3h. 

_ ] 

3t 


(2.19) 


^  i-1 

N  3h . ) 

Add  and  subtract 

i-1 

^  J  -1  J 

E  — i 

.  t'  T  3t 
3=1+1  3 

to  obtain 
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(2.20) 


3A.  8h. 

1  .  1 

3t  '  at 


'  i 

(  N  ah . ) 

'  i-1 

(  N  ah . ) 

j=i  ^  V 

r  _ 1 

.  r  .at 

V  3=1+1  j 

^  J--L  y 

L  O 

.  .  at 

3=1  7 

The  purpose  of  the  preceding  mathematical  manipulations  is  now  apparent. 

Both  the  kinetic  energy  and  available  potential  energy  relations  have  been 

written  in  terms  h.  and  time  derivatives  of  h..  Furthermore,  the  various 

i  i  ^ 

terms  have  been  grouped  to  show  the  transfers  of  energy  between  layers.  This 

grouping  of  terms  for  a  given  hydrodynamic  model  is  usually  neither  obvious 

nor  unique. 


The  energy  transformations  are  more  apparent  when  the  relations  for 
kinetic  energy,  (2.16),  and  available  potential  energy,  (2.20),  are  written 
for  the  layers  above,  i-1,  and  below,  i+1,  the  layer  i  as  follows: 


aK.  , 

1-1 

at 


ah. 


♦  ‘'i-iPi-i)  -  Pi-1  ^ 

(1)  (2) 


(2.21) 


aA.  , 

1-1 


at  ■  Pi-1  at 
(2) 


ah.  , 

1-1 


i-1 

z:  p.gh 
U=i  ^  ^ 


\ 

^  N 

ah.) 

/ 

E 

3 

^  j  =  i 

3t 

/ 

> 

\ 

i-2 

I  p.gh. 

U  =  1  ^  ^ 


N  ah .  1 

j;  _ 12 

.  .  .at 
J=l-1  J 


(3) 


(4) 


(2.22) 


aK, 

^t 


ah. 


-  =  -  ^"v.  [k.  +  h.p.l  -  p. 

1  (  1  i^i)  *^1  at 

(5)  (6) 


(2.23) 


aA. 


ah. 

1 


at  "  Pi  at 


(6) 


r  p.gh 

j=i  ^  ^ 


N 
Z 

I  j=i+l 

(7) 


ah-.^ 

3 

at 


i-1 

I  p.gh. 

U=i  ^  ^ 


N  ah.) 
V  _ i 


at 


(3) 


(2.24) 


aK.  , 

1+1 

at 


-  ^-v.  .  (k.  .  +  h.  .p.  .  I 

1+1  (  1+1  i+i^^i+i; 

(8) 


ah.  , 

1+1 


Pi+i  at 

(9) 


(2.25) 


aA.  , 

1  +  1 

at 


=  p 


ah.  , 

1+1 


i+1  at 

(9) 


i+1 

Z  p.gh. 
V  j=l  ^  ^ 


N  ah.) 


j  =  i  +  2 


at 


(10) 
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(2.26) 


'  i  ^ 

'  N  ah . " 

i-1 

V  3 -J-  / 

.  ^  IT^ 
j=i+i  J 

(7) 


Consider  the  energy  transformations  for  layer  i.  Term  6  is  the  transfer 
between  kinetic  energy  and  available  potential  energy  within  the  layer.  Term 
7  in  (2.24)  and  (2.26)  is  the  transfer  of  available  potential  energy  between 
layer  i  and  layer  i+1.  Term  3  in  (2.22)  and  (2.24)  is  the  transfer  of 
available  potential  energy  between  layer  i-1  and  layer  i.  It  follows  that 
term  4  in  (2.22)  is  the  transfer  of  available  potential  energy  between  layer 
i-1  and  layer  i-2  and  that  term  10  in  (2.26)  is  the  transfer  of  available 
potential  energy  between  layer  i+1  and  layer  i+2.  For  completeness,  it  should 
be  noted  that  term  2  in  (2.21)  and  (2.22)  and  term  9  in  (2.25)  and  (2.26) 
represent  the  transfer  between  kinetic  and  available  potential  energy  in 
layers  i-1  and  i+1,  respectively.  Finally,  terms  1,  5,  and  8  are  the  energy 
flux  divergence  terms  for  layers  i-1,  i,  and  i+1,  respectively.  These 
relations  are  shown  in  the  block  diagram  in  Figure  2.  Note  that  only  layer  i 
has  all  transfers  marked. 

The  preceding  discussion  indicates  that  the  kinetic  energy  and  available 
potential  energy  relations,  (2.16)  and  (2.20),  may  be  written  and  interpreted 
as  follows: 


8K. 
_ 1 

at 


4  h.p.) 

(1) 


ah. 
_ 1 

i  at 

(2) 


aA. 

1 

at 


ah. 

^  i  ' 

N  ah . ) 

^  i-1 

r  N  ah . ) 

Pi  at'  * 

r  p.Bh. 

r  _ i 

.  ^  at 

j=i+i  j 

- 

^  3-*^  / 

r  _ 1 

.^.  at 

1  j=i  J 

(2)  (3)  (4) 


(2.27) 


(2.28) 


Term  1  is  the  energy  flux  divergence  for  layer  i.  Term  2  is  the  transfer 
between  kinetic  energy  and  available  potential  energy  in  layer  i.  Term  3  is 
the  transfer  of  available  potential  energy  between  layer  i  and  layer  i+1  while 
term  4  is  the  transfer  between  layer  i  and  layer  i-1. 


The  available  potential  energy  relation,  (2.28),  shows  transfers  of 
energy  with  layers  directly  above  and  below  the  layer  i.  What  happens  when 
layer  i  is  the  top  or  bottom  layer?  For  the  top  layer,  i=l,  (2.28)  takes  the 
form 


aA^  ah^  N  ah 

yr  "  ^1  IF  ^i^^i  at 

3 


(2.29) 
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The  second  term  on  the  RHS  of  (2.29)  is  the  transfer  of  available  potential 
energy  between  layers  1  and  2.  As  there  is  no  layer  above  layer  1,  there  is 
no  term  for  the  transfer  of  available  potential  energy  between  layer  1  and  the 
layer  above  layer  1.  For  the  bottom  layer,  i=N,  (2.28)  takes  the  form 


3t 


=  P 


3h 
N  dt 


N 


N-1 

I  P.gh. 

j=l  ^  J 


at 


(2.30) 


The  second  term  on  the  RHS  of  (2.30)  is  the  transfer  of  available  potential 
energy  between  the  bottom  layer  N  and  the  layer  above,  N-1.  As  there  is  no 
layer  below  the  bottom  layer,  there  is  no  term  for  the  transfer  of  available 
potential  energy  between  layer  N  and  the  layer  below  layer  N. 

The  total  energy  per  unit  area  in  a  layer  is  K^+A^.  Adding  (2.27)  and 
(2.28)  gives  the  local  time  derivative  of  the  total  energy  per  unit  area  in  a 
layer  as 


f  ^  1 

'  N 

ah.  ) 

It  (‘'i  *  *i] 

1  ■  - 

(k.  *  h.p.]  . 

Z 

. j=i+l 

3 

at 

> 

(1)  (3) 


i-1 

I  p.gh. 

j=l  ^  ^ 


N  ah.  ) 


J=1 


at 


(4) 


(2.31) 


Here,  term  1  is  the  energy  flux  divergence  while  terms  3  and  4  represent  the 
exchange  of  total  energy  per  unit  area  between  the  layer  i  and  the  layers 
below  and  above,  respectively. 

If  the  total  energy  per  unit  area  in  a  layer,  K^+A^,  is  summed  over  all 
layers,  the  sum  is  called  the  total  energy  per  unit  area  of  the  system. 
Summing  (2.31)  over  all  N  layers  gives  the  local  time  derivative  of  the  total 
energy  per  unit  area  for  the  system  as 


N 

Z 

i=l 


at 


K 


+  A. 
1 


*iPi) 


(2.32) 


Only  flux  divergence  terms  appear  on  the  RHS  of  (2.32);  there  is  no  longer  any 
transfer  of  energy  between  layers.  Nov  integrate  (2.32)  over  the  horizontal 
area  of  the  system  to  get  the  local  time  derivative  of  the  total  energy  in  the 
entire  system: 


N 

E 

i=l 


at 


(k.  .  aJ 


dxdy 


Z  ^-v.  (k.  +  p.h.l  dxdy 
.  .  ill  11) 

1=1 
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N 

Z 

i=l  . 


S. 

1 


K./h. 
I  1  1 


-> 

V. 

1 


•n. 

1 


dS. 

1 


(2.33) 


where  S.  is  the  total  surface  of  the  side  walls  in  layer  i,  is  the  unit 
vector  normal  to  the  side  walls  in  layer  i,  and  dS.  represents  the 
differential  area  of  the  side  walls  in  layer  i.  If  the  side  walls  are 
considered  to  be  solid,  there  can  be  no  flux  of  energy  through  them  so  (2.33) 
becomes 


Z  +  A^j  dxdy  =  0  .  (2.34) 


This  relation  merely  states  that  the  total  energy  is  conserved  for  an 
inviscid  multilayer  hydrodynamic  model  with  no  external  sources  or  sinks. 

The  energy  specifications  and  relations  in  the  layer-layer  formulation  for  the 
finite  depth  model  are  given  for  the  one-,  two-,  and  three-layer  cases  in 
sections  C.8,  C.9,  and  C.IO,  respectively. 


It  K  "  ''i)  ■  It 


For  some  studies  the  total  kinetic  energy  and  the  total  available  potential 
energy  are  treated  separately.  The  local  time  derivative  of  the  total  kinetic 
energy  for  the  system  is  found  by  summing  (2.27)  over  all  N  layers  and  then 
integrating  over  the  horizontal  area: 
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or,  for  solid  side  walls. 
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(2.35) 


(2.36) 


Similarly  the  local  time  derivative  of  the  total  available  potential  energy 
for  the  system  is  found  by  summing  (2.28)  over  all  N  layers  and  then 
integrating  over  the  horizontal  area: 


N  aA. 
i=l 


N  ah. 

2:  p.  ^  dxdy 


(2.37) 
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If  (2.36)  and  (2.37)  are  added,  the  result  is  (2.34)  which  shows  that  energy 
is  transferred  between  kinetic  and  available  potential  forms  even  though  the 
total  energy  for  the  basin  is  conserved.  See  Semtner  and  Mintz  (1977)  and 
Robinson  et  al.  (1979)  for  further  examples  of  basin-integrated  energetics. 
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Part  B,  Section  3 


In  sections  B.l  and  B.2  the  available  potential  energy  per  layer  has  been 
used  in  determining  energy  transfers.  However,  there  is  another  way  to  treat 
the  available  potential  energy  that  will  be  derived  in  this  section. 

From  (1.22),  the  available  potential  energy  per  unit  area,  A^,  in  layer  i 
can  be  written  as 
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(3.1) 


The  total  available  potential  energy  per  unit  area  for  the  entire  system  is 
found  by  summing  (3.1)  over  the  N  layers  to  get 
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(3.3) 


Examination  of  (3.3)  reveals  that  available  potential  energy  is _ associated 
with  each  interface  deviation,  \]^.  Define  the  available  potential  energy 
associated  with  interface  deviation,  such  that 
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If  (3.4)  is  summed  over  N  layers,  the  result  is 
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(3.6) 
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Comparison  of  (3.3)  and  (3.6)  shows  that  the  total  available  potential  energy 
per  unit  area  in  the  system  is  the  same  no  matter  that  formulation  of 
available  potential  energy  is  used:  layer,  A^,  or  interface,  A^. 

Take  8/9t  of  (3.4)  and  (3.5)  to  get 
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(3.8) 

From  (2.16),  the  change  w.r.t.  time  of  kinetic  energy  per  unit  area  in  a  layer 
is 
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where,  from  (2.3)  and  (2.5), 
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111  1+1 
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(3.10) 


(3.11) 


The  interface  available  potential  energy.  A.,  is  given  in  (3.4)  -  (3.8)  in 
terms  of  the  interface  deviation,  t).  ,  so  use  (3.10)  in  the  last  term  on  the 
RHS  of  (3.9)  to  eliminate  h^^.  The  result  is 
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Now  evaluate  p^  -  by  using  (3.11)  twice  to  get 

Pi  -  Pi_i  =  Pig  hi  (pi  -  pjg  hj  -  Pi_ig  hi  (Pi_i  -  pjg  h.  (3.13) 
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or 
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(3.14) 


Sum  (3.10)  from  layer  1  to  layer  i-l  to  get 
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Use  (3.16)  in  (3.14)  to  find 
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(3.18) 


Use  (3.18)  in  (3.7)  to  obtain 
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Both  kinetic  energy  and  interface  available  potential  energy  relations 
have  been  written  in  terms  of  the  interface  deviation,  and  its  derivative 

w.r.t.  time.  It  is  useful  to  consider  the  energy  relations  (3.19)  and  (3.12) 
for  adjacent  interfaces  and  layers.  For  i-l,  i,  and  i+1,  (3.19)  and  (3.12) 
become 
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(9) 

(8) 

(10) 

The  transfer  between  kinetic  energy  in  a  given  layer  and  the  available 
potential  energy  for  the  interface  above  the  given  layer  is  given  by  terms  1, 
5,  and  8  for  layers  i-1,  i,  and  i+1,  respectively.  Similarly,  the  transfer 
between  kinetic  energy  in  a  given  layer  and  the  available  potential  energy  for 
the  interface  below  the  given  layer  is  given  by  terms  2,  4,  7,  and  10  for 
layers  i-2,  i-1,  i,  and  i+1,  respectively.  Terms  3,  6,  and  9  are  the  energy 
flux  divergences  for  layers  i-1,  i,  and  i+1,  respectively.  These  relations 
are  shown  in  the  block  diagram  in  Figure  3. 

The  preceding  discussion  indicates  that  kinetic  energy  and  interface 
available  potential  energy  relations  may  be  written  and  interpreted  as 
follows: 
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Terms  2  and  3  are  the  transfers  between  the  kinetic  energy  in  layer  i  and  the 
available  potential  energy  for  the  interfaces  above  and  below  layer  i, 
respectively.  Term  4  is  the  transfer  between  the  available  potential  energy 
for  interface  i  and  the  kinetic  energy  in  layer  i-1,  the  layer  above  interface 
i.  Term  1  is  the  usual  energy  flux  divergence  associated  with  layer  i.  What 
happens  to  these  relations  for  the  top  and  bottom  layers?  For  the  top  layer, 
(3.27)  takes  the  form 


8t“  "  IF' 


i=l  . 


(3.28) 


As  there  is  no  layer  above  interface  1,  no  term  exists  for  the  transfer  of 
energy  to  that  layer.  For  the  bottom  layer,  i=N,  (3.26)  becomes 


at 


Pn  at 


i=N  . 


(3.29) 


As  there  is  no  interface  below  layer  N,  there  is  no  term  for  the  transfer  of 
any  energy  associated  with  such  an  interface. 

The  local  time  derivative  of  the  total  kinetic  energy  per  unit  area  is 
found  by  summing  (3.26)  over  N-1  layers  and  adding  (3.29)  to  get 
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(3.30) 


The  local  time  derivative  of  the  total  available  potential  energy  per 
unit  area  is  found  by  summing  (3.27)  from  i=2  to  i=N  and  adding  (3.28)  to  find 
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or,  rearranging  terms. 
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Adding  (3.30)  and  (3.32)  gives  the  local  time  derivative  of  the  total  energy 
per  unit  area  as 
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(3.33) 
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1^1 


Compare  this  energy  relation  to  (2.32)  that  is  the  change  in  the  total  energy 
per  unit  area  for  the  system  in  the  layer-layer  formulation.  In  both 
formulations  only  energy  flux  divergence  terms  contribute  to  a  change  w.r.t. 
time  in  the  total  energy  per  unit  area. 

If  (3.33)  is  integrated  over  the  horizontal  area  of  the  basin,  the  result 
is 
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(3.34) 


for  solid  side  walls  as  found  in  (2.34)  for  the  layer-layer  formulation  of 
energetics . 


The  energy  specifications  and  relations  in  the  layer-interface 
formulation  for  the  finite  depth  model  are  given  in  sections  D.13,  D.14,  and 
D.15  for  the  one-,  two-,  and  three-layer  cases,  respectively. 
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Part  B,  Section  4 


The  major  difference  between  the  layer-layer  and  the  layer-interface 
formulations  in  sections  B.2  and  B.3  lies  in  the  grouping  of  terms.  The 
total  kinetic  energy,  summed  over  all  layers,  remains  the  same.  The  total 
available  potential  energy,  whether  summed  over  all  layers  in  the  layer 
formulation  or  summed  over  all  interfaces  in  the  interface  formulation, 
remains  the  same  for  the  system  as  shown  in  (3.3)  and  (3.6).  The  question 
naturally  arises  as  to  whether  there  are  other  ways  of  formulating  the  energy 
transfers.  In  this  section,  a  third  way  of  formulating  energy  transfers  for 
the  finite  depth  hydrodynamic  model  is  derived.  It  is  included  here  as  an 
example  of  an  energy  formulation  that  is  correct  mathematically,  but 
unsatisfactory  in  terms  of  the  physical  mechanisms  of  transfer  represented. 

Consider  the  layer-level  formulation  in  section  B.3.  The  relevant  energy 
transfer  relations,  (3.26)  and  (3.27),  are 
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Now  add  and  subtract  the  term  p^_j  to  (4.1)  to  derive 
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It  is  useful  to  consider  (4.2)  and  (4.3)  for  adjacent  interfaces  and  levels. 
For  i-1,  i,  and  i+1,  these  relations  become 
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Terms  1,  5,  and  8  represent  the  transfer  between  the  kinetic  energy  in  a  layer 
and  the  available  potential  energy  associated  with  the  interface  directly 
above  that  layer  for  layers  i-l,  i,  and  i+1,  respectively.  Term  3  is  the 
transfer  of  kinetic  energy  between  layers  i-l  and  i  and  term  7  is  the  same 
transfer  between  layers  i  and  i+1.  Terms  4  and  10  are  the  transfer  of  kinetic 
energy  between  layers  i-l  and  i-2  and  layers  i+1  and  i+2,  respectively.  Terms 
2,  6,  and  9  represent  the  energy  flux  divergence  per  unit  area  in  layers  i-l, 
i,  and  i+1,  respectively.  These  relations  are  shown  in  the  block  diagram  in 
Figure  4. 

The  preceding  discussion  indicates  that  the  kinetic  energy  and  available 
potential  energy  relations,  (4.2)  and  (4.3)  can  be  written  and  interpreted  as 
follows ; 
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Term  1  is  the  energy  flux  divergence  for  layer  i.  Term  2  is  the  transfer 
between  kinetic  energy  in  layer  i  and  the  available  potential  energy 
associated  with  the  interface  i  above  layer  i.  Term  3  is  the  exchange  of 
kinetic  energy  between  layer  i  and  the  layer  above,  i-1.  Term  4  is  the 
exchange  of  kinetic  energy  between  layer  i  and  the  layer  below,  i+1. 

The  formulation  of  energy  transfers  in  this  section  is  unsatisfactory  as 
the  physical  mechanisms  depicted  are  unrealistic.  As  shown  in  Figure  4, 
kinetic  energy  in  a  layer  is  transferred  to  and  from  available  potential 
energy  for  the  interface  above  the  layer.  There  is  no  reason  such  an  exchange 
cannot  take  place  between  a  given  layer  and  the  interface  below,  as  shown  in 
Figure  3,  yet  no  mechanism  for  such  a  transfer  occurs  in  the  formulation  for 
this  section.  Furthermore,  this  formulation  allows  for  the  transfer  of^ 
kinetic  energy  between  layers.  However,  the  inviscid  ocean  model  described  in 
the  Appendix  does  not  permit  the  direct  transfer  of  mass  from  one  layer  to 
another,  so  there  is  no  physical  mechanism  to  allow  the  direct  exchange  of 
kinetic  energy  between  layers.  '  These  two  shortcomings  indicate  that  the 
formulation  in  this  section,  although  mathematically  correct,  is  not 
satisfactory  in  terms  of  physical  mechanisms  needed  for  this  ocean  model. 
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Part  B,  Section  5 


In  this  section  the  specific  equations  for  energy  transfer  between  layers 
and  within  a  layer  are  found  in  a  layer-layer  formulation  for  a  multilayer 
reduced-gravity  model.  The  reduced-gravity  model  contains  one  more  assumption 
than  is  given  in  the  Appendix:  the  bottom  layer,  N,  is  assumed  to  be  of 
infinite  depth.  How  does  this  additional  assumption  alter  the  energy 
relations  found  in  section  B.2? 

From  section  B.l,  (1.1),  the  kinetic  energy  per  unit  area  in  layer  N,  Kj^,  is 
written  as 
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where,  from  (2.4) 

Substitute  (5.2)  in  (5.1)  to  derive 
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Now  as  H  “,  K^->  ®,  unless  u  and  v^->  0.  It  is  not  physically  reasonable  to 
have  infinite  kinetic  energy  so  the  velocity  in  the  bottom  layer  is  set  to^ 
zero  in  a  reduced-gravity  model.  It  follows  that  0  and  remains  zero,  i.e. 
8K^/9t  =  0.  Hence,  the  bottom  layer  is  said  to  be  inactive  in  a  reduced- 
gravity  model. 

The  starting  equations,  derived  in  the  Appendix  and  listed  in  Table  1, 

are 
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These  equations  for  layers  i<N  do  not  depend  on  the  depth  of  layer  N  and  so 
remain  valid  for  layers  i<N  in  the  reduced-gravity  model.  What  happens  to 
these  equations  for  the  bottom  layer,  i=N,  in  the  reduced-gravity  model? 

Consider  first  the  pressure  equation,  (5.7).  From  (A. 7),  the  pressure, 
Pj,  can  be  written  as 

p.(x,y,z,t)  =  p.(x,y,t)  -  p.gz.  (5.8) 

If  (5.8)  is  used  in  (5.7)  to  eliminate  p^,  the  result  for  the  bottom  layer  is 
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so  (5.9)  may  be  written  as 
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(5.10) 


Thus,  Pj^  does  not  depend  on  the  depth  of  the  bottom  layer  and  (5.10)  may  be 
used  for  the  reduced-gravity  model. 

Consider  next  the  momentum  equations  (5.4)  -  (5.5)  for  the  bottom  layer 
in  the  reduced-gravity  model.  The  velocity  is  zero  in  the  bottom  layer  so 
these  equations  become 
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Thus,  for  the  reduced-gravity  model,  the  pressure  in  the  bottom  layer  is  not  a 
function  of  x  and  y,  so  p„  is  a  function  of  time  only.  Next,  average  (5.10) 
over  the  area  of  the  model: 
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(5.13) 
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where  the  angle  brackets  indicate  an  average  over  the  basin  area.  From  (2.3) 
the  relation  between  and  h^  is  given  as 


h.  =  t). 
1  1 


n.  ,  +  H. 
'i+l  1 


(5.14) 


The  basin  area-average  of  (5.14)  is 


=H. 


i<N 


(5.15) 


because  H.  is  defined  as  the  basin  area-averaged  thickness  of  layer  i. 
Furthermore,  h-  is  the  difference  between  the  actual  and  basin  area-averaged 
position  of  the  top  of  layer  i.  It  follows  that 

[.  =  0  l<i^  •  (5.16) 


Now  the  mass  in  layer  i,  is  given  by 

M.  =  p.  •^h.  \  X  Area  =  p.H.  x  Area  i<N  .  (5.17) 

For  the  multilayer  hydrodynamic  models  considered  here,  mass  in  a  layer  is 
assumed  to  be  conserved  so  M.  and  therefore  do  not  vary  with  time.  If 
(5.15)  and  (5.16)  are  used  in  (5.13),  the  result  is 
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(5.18) 


which  is  a  constant  for  the  system  and  not  a  function  of  x,  y,  or  t.  Note 
that  (5.18)  is  correct  for  both  the  finite  depth  and  the  reduced-gravity 
models.  However,  for  the  reduced-gravity  model,  it  follows  from  (5.11)  and 
(5.12)  that 


Pn  -  “Ipn^ 

and  from  (5.10),  (5.18),  and  (5.19), 
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(5.20) 


Thus,  in  the  reduced-gravity  model  Pj^  is  also  not  a  function  of  time.  Sum 
(5.14)  over  N-1  layers  to  obtain 
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Use  (5.21)  to  eliminate  in  (5.20);  the  result  is 


N-l  , 
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Take  3/3t  of  (5.22)  to  find  that 
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The  deriviative  v.r.t.  time  of  (5.2) 
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is 


(5.21) 


(5.22) 


(5.23) 


3^N 
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(5.24) 


that  is  valid  even  if  ->  So  (5.23)  can  be  written  as 


3n^  N-l 


Note  that  (5.6)  has  been  obtained  in  the  Appendix  by  integrating  over  the 
layer  depth.  Thus,  in  a  reduced-gravity  model,  (5.6)  does  not  apply  to  the 
bottom  layer,  but  (5.24)  or  (5.25)  is  used  to  evaluate  the  change  v.r.t.  time 
of  the  interface  deviation, 

The  preceding  discussion  reveals  that  the  assumption  of  infinite  depth 
for  the  bottom  layer  in  a  multilayer  hydrodynamic  model  has  three  direct 
consequences.  First,  the  velocity  in  the  bottom  layer  is  zero.  Second,  the 
pressure  in  bottom  layer  is  independent  of  x,  y,  and  t.  Third,  as  shown  in 
(5.24),  the  interface  deviation  at  the  top  of  the  bottom  layer,  hj.,  cannot  be 
specified  independently,  but  is  directly  coupled  to  the  sum  of  all  other 
interface  deviations. 

Consider  now  the  energy  relations  for  a  reduced-gravity,  hydrodynamic 
model.  For  layers  1  to  N-l,  (5.4)  -  (5.7)  are  the  same  equations  used  to 
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3t 
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(5.25) 
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derive  the  energy  relations  in  section  B.2  for  the  layer-layer  formulation  for 
the  finite  depth  model.  These  energy  relations  are,  from  (2.27)  and  (2.28), 


— i  =  -^‘v.fK.  +h.p.l-  p. 
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For  the  bottom  layer,  the  discussion  following  (5.3)  has  shown  that 


at 


=  0  . 


(5.28) 


The  general  expression  for  the  local  time  derivative  of  available  potential 
energy  is  derived  in  section  B.l  and  does  not  depend  on  the  layer  depth,  only 
on  interface  deviation.  Thus,  this  relation,  (1-24),  is  valid  for  the 
reduced-gravity  model  and  is,  for  layer  N, 
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From  (5.22),  can  be  written  as 
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(5.30) 


If  (5.14)  is  used  in  (5.30),  the  result  is 
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(5.31) 


Therefore,  in  the  reduced-gravity  model,  the  interface  deviation  at  the  top  of 
the  bottom  layer  is  determined  completely  by  a  sum  of  the  other  interface 
deviations  multiplied  by  density  differences.  This  result,  (5.31),  is  a 
direct  consequence  of  the  fact  that  the  pressure  in  the  bottom  layer  is 
independent  of  horizontal  position  and  time,  and  it  does  not  hold  for  the 
finite  depth  model. 
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Use  (5.30)  and  (5.24)  in  (5.29)  to  eliminate  the  result  is 
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With  the  aid  of  (5.22),  (5.32)  can  be  written  as 
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A  mathematically  consistent  set  of  energy  relations  has  been  derived  for  a 
reduced-gravity  model:  (5.26)  and  (5.27)  for  layers  1  to  N-l  and  (5.28)  and 
(5.33)  for  layer  N.  However,  there  is  a  problem  with  the  interpretation  of 
the  energy  terms.  To  illustrate  this  problem,  note  that  for  i=N-l,  (5.27) 
becomes 
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For  the  layer-layer  formulation,  available  potential  energy  transfer  occurs 
between  adjacent  layers  and  transfer  between  kinetic  energy  and  available 
potential  energy  occurs  within  a  layer.  In  (5.34)  the  first  term  on  the  RHS 
is  transfer  between  kinetic  energy  and  available  potential  energy  while  the 
second  term  is  the  transfer  of  available  potential  energy  between  layers  N  and 
N-l.  In  the  bottom  layer  there  is  no  transfer  between  kinetic  energy  and 
available  potential  energy;  both  terms  on  the  RHS  of  (5.33)  correspond  to  a 
transfer  of  available  potential  energy  between  layers  N  and  N-l.  It  is 
obvious  that  the  terms  for  the  transfer  of  available  potential  energy  between 
layers  N  and  N-l  in  (5.33)  and  (5.34)  do  not  match  as  they  should.  It  is 
necessary  to  regroup  the  factors  in  (5.34)  and,  therefore,  (5.27). 

Define  a  new  quantity,  p^,  as 
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(5.36) 
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or 
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For  i=N,  (5.37)  becomes 
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which  is  the  same  as  (5.33).  For  i=N-l,  (5.37)  becomes 
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Nov,  terms  1  and  2  in  (5.39)  match  those  on  the  RHS  of  (5.38),  so  these  terms 
are  the  transfer  of  available  potential  energy  between  layers  N  and  N-1. 

Terms  4  and  5  of  (5.39)  give  the  transfer  of  available  potential  energy 
between  layers  N-1  and  N-2.  Term  3  should  correspond  to  the  transfer  between 
kinetic  energy  and  available  potential  energy  in  layer  N-1.  Does  such  a 
matching  term  exist  in  the  kinetic  energy  relation.  Use  (5.35)  in  (5.26)  to 
replace  p^: 
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or,  for  layer  N-1, 
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(5.41) 
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Term  3  appears  in  (5.39)  and  (5.41)  and  is  the  transfer  between  kinetic  and 
available  potential  energy  within  the  layer  N-1.  For  N-2,  (5.37)  and  (5.40) 
become 
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These  relations  are  shown  in  the  block  diagram  in  Figure  5. 

The  preceding  discussion  indicates  that  kinetic  energy  and  available  potential 
energy  relations  for  the  reduced-gravity  model  can  be  written  and  interpreted 
as  follows: 
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Term  1  is  the  modified  energy  flux  divergence.  Term  2  is  the  transfer  between 
kinetic  energy  and  available  potential  energy  within  a  layer.  Term  3  is  the 
transfer  of  available  potential  energy  between  layer  i  and  the  layer  below. 
Term  4  is  the  transfer  of  available  potential  energy  between  layer  i  and  the 
layer  above. 

Note  that  v.,  and  p„  are  zero,  so  for  i=N,  (5.44)  becomes  (5.36)  or 
N  N 


However,  it  is  desirable  to  apply  (5.44)  only  for  layers  i<N  because  (5.6)  has 
been  used  to  get  (5.44)  and  (5.6)  is  valid  only  for  i<N. 

The  local  time  derivative  of  the  total  kinetic  energy  per  unit  area  is 
found  by  summing  (5.44)  over  N-1  layers  and  adding  (5.45)  to  get 
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Summing  (5.46)  over  N  layers  shows  that  the  local  time  derivative  of  the  total 
available  potential  energy  per  unit  area  is 
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Adding  (5.48)  and  (5.49)  gives  the  local  time  derivative  of  the  total  energy 
per  unit  area  as 
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A  comparison  of  (5.50)  and  (2.32),  that  is  the  change  w.r.t.  time  of  the  total 
energy  per  unit  area  of  the  system  for  finite  depth  model  shows  that  the  RHS 
of  both  equations  appear  to  contain  energy  flux  divergence  terms.  Note, 
however,  that  (5.48)  is  summed  over  only  N-1  layers;  the  bottom  layer  does  not 
contribute  an  energy  flux  divergence  term  in  the  reduced-gravity  model. 
Furthermore,  energy  flux  divergence  terms  for  the  reduced-gravity  model  use 
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p. ,  not  p.,  as  for  the  finite  depth  model.  Use  (5.35)  to  eliminate  p^  in 
(i.48).  i'he  result  is 
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It  is  apparent  from  (5.52)  that  the  energy  flux  divergence  term  for  the  bottom 
layer  that  is  present  in  the  finite  depth  model,  (2.32),  has  been  replaced  by 
p  /3t  in  the  reduced-gravity  model.  Now  integrate  (5.52)  over  the 

horizontal  area  of  the  system  to  get  the  local  time  derivative  of  the  total 
energy  for  the  entire  system  of  the  reduced-gravity  model: 
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where  S.  is  the  total  surface  of  the  side  walls  for  layer  i,  n.  is  the  unit 
vector  normal  to  the  side  walls  and  dS. represents  the  differential  area  of  the 
side  walls  for  layer  i.  Now  from  (5.1b),  ~  that  last  term  on  the 

RHS  of  (5.54)  is  zero.  If  the  side  walls  are  solid,  there  is  no  flux  through 
them  and  the  first  term  on  the  RHS  of  (5.54)  vanishes,  so 


E  +  A^jdxdy  =  0  .  (5.55) 

This  relation  merely  states  that  the  total  energy  is  conserved  for  an 
inviscid  multilayer  reduced-gravity  hydrodynamic  model  with  no  external 
sources  or  sinks.  Thus,  as  a  comparison  of  (2.34)  and  (5.55)  shows,  the 
additional  assumption  of  an  infinitely  deep  bottom  layer  does  not  change  the 
fact  that  total  energy  is  conserved  for  a  system. 


at 
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The  energy  specifications  and  relations  in  the  layer-layer  formulation 
for  the  reduced-gravity  model  are  given  in  sections  C.ll  and  C.12  for  the  one 
and  two  active  layer  cases,  respectively. 
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Part  B,  Section  6 


In  this  section  the  specific  equations  for  transfer  of  kinetic  energy  and 
available  potential  energy  are  found  in  a  layer-interface  formulation  for  an 
inviscid  multilayer  reduced-gravity  hydrodynamic  model.  The  reduced-gravity 
model  is  discussed  extensively  in  section  B.5  and  the  layer-interface 
formulation  is  treated  in  section  B.3.  Both  sections  should  be  perused  before 
reading  further  in  this  section. 

The  starting  equations,  derived  in  the  Appendix,  remain  valid  for  the 
first  N-1  layers.  Therefore,  the  same  energy  relations  are  valid  for  these 
layers.  These  relations,  (3.26)  -  (3.28),  are 
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In  a  reduced-gravity  model  both  the  velocity,  Vj^,.and  the  kinetic  energy, 
are  zero  in  the  bottom  layer  while  the  pressure,  p^^,  is  constant.  From 
(5.22),  the  pressure  can  be  written  as 
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Section  B.5  shows  that  Bhj^  /at  cannot  be  determined  from  (A. 17),  which  is  not 
valid  in  the  bottom  layer,  but  it  can  be  found  from  (5.25)  as 
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The  available  potential  energy  relation,  (6.2),  was  derived  in  section  B.3  and 
does  not  depend  on  the  depth  of  the  bottom  layer.  So,  for  ir:N,  (6.2)  becomes 
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Because  the  kinetic  energy  is  constant  and  zero  in  the  bottom  layer,  the 
kinetic  energy  relation  for  the  bottom  layer  is 


(6.8) 


A  mathematically  consistent  set  of  energy  relations  has  been  derived  for 
a  reduced-gravity  model  in  the  layer-interface  formulation:  (6.1)  -  (6.3)  for 
layers  1  to  N-1  and  (6.7)  and  (6.8)  for  layer  N.  However,  as  is  the  case  for 
the  layer-layer  formulation,  there  is  a  problem  in  interpreting  the  energy 
transfer  terms.  To  illustrate  this  problem,  note  that  for  i=N-l,  (6.1)  and 
(6.2)  become 
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Term  2  appears  to  be  the  transfer  between  available  potential  energy  in 
interface  N  and  the  kinetic  energy  in  layer  N-1.  Term  4  appears  to  be  the 
transfer  between  the  kinetic  energy  in  layer  N-1  and  the  available  potential 
energy  in  interface  N-1.  What  is  term  1  in  (6.7)?  It  must  be  a  term  to 
transfer  energy  upward,  but  there  is  no  matching  term  in  (6.9)  and  (6.10).  It 
is  obvious  that  the  terms  need  to  be  regrouped  in  (6.1),  (6.2)  and  (6.7). 

In  section  B.5,  a  new  quantity,  p^,  is  defined  in  (5.35)  as 
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where  Pjj  is  a  constant  in  the  reduced-gravity  model  and  =  0. 
Substituting  (6.11)  into  (6.1)  and  (6.2)  gives 
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For  i=N“2,  (6.14)  and  (6.12)  become 
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For  i=N-l,  (6.14)  and  (6.12)  become 
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(6.15) 

(6.16) 
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Substituting  (6.11)  into  (6.7)  and  (6.8)  yields 
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Terms  5  and  9  are  the  modified  energy  flux  divergences  for  layers  N-2  and  N-1, 
respectively.  Terms  1  and  7  are  the  transfer  between  kinetic  energy  in  a 
given  layer  and  available  potential  energy  of  the  interface  directly  above  for 
layers  N-2  and  N-1,  respectively.  Terms  6  and  10  give  the  transfer  between 
kinetic  energy  in  a  given  layer  and  available  potential  energy  of  the 
interface  directly  below  for  layers  N-2  and  N-1,  respectively.  Terms  8  and  4 
represent  the  transfer  of  available  potential  energy  between  a  given  interface 
and  the  interface  above  for  interfaces  N  and  N-1,  respectively.  These 
relations  are  shown  in  the  block  diagram  in  Figure  6. 

A  comparison  of  Figure  6  to  Figure  3  shows  that  the  introduction  of  the 
reduced-gravity  assumption  has  changed  the  types  of  energy  transfer  possible 
in  the  layer-interface  formulation.  In  Figure  3,  there  is  no  transfer  of 
available  potential  energy  between  interfaces.  It  seems  physically  reasonable 
for  energy  to  be  transferred  between  a  layer  and  an  adjacent  interface.  The 
interface  deviations  are  not  directly  related  for  the  model  in  section  B.3. 
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In  Figure  6,  there  is  a  transfer  of  available  potential  energy  between 
consecutive  interfaces.  This  coupling  occurs  in  the  reduced-gravity  model 
because  the  interface  deviations  are  directly  related  through  (6.6)  which  is 
valid  only  for  the  reduced-gravity  model. 

The  preceding  discussion  indicates  that  kinetic  energy  and  available 
potential  energy  relations  for  the  reduced-gravity  model  in  layer-interface 
formulation  can  be  written  and  interpreted  as  follows; 
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(6.21) 


(6.22) 


(6.23) 


Term  2  is  the  transfer  between  the  kinetic  energy  in  layer  i  and  the  available 
potential  energy  for  the  interface  i,  directly  above.  Term  3  and  term  4  are 
the  same  type  of  transfer  for  layer  i  and  interface  i+l  and  for  layer  i-1  and 
interface  i,  respectively.  Terms  5  and  6  represent  a  transfer  of  available 
potential  energy  from  interface  i  to  the  interface  above  and  below, 
respectively.  Finally,  terra  1  is  the  modified  energy  flux  divergence. 

Summing  (6.21)  over  N-1  layers  and  adding  (6.22)  gives  the  time 
derivative  of  the  total  kinetic  energy  per  unit  area  for  the  reduced-gravity 
model  as 
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(6.24) 


Summing  (6.23)  over  N  interfaces  gives  the  time  derivative  of  the  total 
interface  available  potential  energy  per  unit  area  for  the  reduced-gravity 
model  as 
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(6.25) 
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or,  rearranging  terms, 
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(6.26) 


Adding  (6.24)  and  (6.26)  yields  the  time  derivative  of  the  total  energy  per 
unit  area  for  the  reduced-gravity  model  as 


i.l 
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(6.27) 

VI  1  / 

in  the 
(5.6), 

layer-interface  formulation.  Eliminate  p^  in  (6.27) 
and  (5.25);  the  result  is 

by  using  (6.11) , 
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(6.28) 

As  discussed  in  section  B.5  following  (5.52)  the  energy  flux  divergence  term 
for  the  bottom  layer  in  the  finite  depth  model,  (2.32),  is  replaced  by  the 
last  term  on  the  RHS  of  (6.28).  Nov  integrate  (6.28)  over  the  horizontal  area 
of  the  basin  to  get  the  local  time  derivative  of  the  total  energy  in  the  basin 
for  the  reduced-gravity  model: 


at 
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(6.29) 


(6.30) 


If  the  side  walls  are  solid,  then  there  is  no  flux  through  the  side  walls  and 
(6.30)  becomes 
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(6.31) 


=  0  . 


This  relation  indicates  that  in  the  layer-interface  formulation  the  total 
energy  in  the  basin  is  conserved  for  the  reduced-gravity  model  as  expected 
from  (5.55). 

The  energy  specifications  and  relations  in  the  layer-interface 
formulation  for  the  reduced-gravity  model  are  given  in  sections  D.16  and  D.17 
for  the  one  and  tvo  active  layer  cases,  respectively. 
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Part  B,  Section  7 


The  energy  relations  developed  in  sections  B.l  -  B.3  for  the  finite  depth 
model  and  in  sections  B.5  -  B.6  for  the  reduced-gravity  model  illustrate  two 
problems  common  to  all  energy  analysis  in  geophysical  fluid  dynamics:  first, 
specification  of  the  energies  present  and,  second,  grouping  of  energy  transfer 
terms  so  that  realistic  physical  mechanisms  are  presented.  Unfortunately,  no 
foolproof  method  exists  for  handling  these  two  problems.  In  fact,  the  work  in 
section  B.4  is  included  to  show  how  easy  it  is  to  develop  a  mathematically 
correct,,  but  physically  unsatisfactory,  set  of  energy  relations. 

Specification  of  the  energies  may  seem  straightforward  at  first, 
particularly  for  hydrodynamic  models  that  use  only  kinetic  and  potential 
energies.  However,  only  that  portion  of  the  potential  energy  that  can  be 
converted  into  kinetic  energy  is  of  interest  in  an  energy  analysis.  Thus,  in 
section  B.l  this  available  potential  energy  is  specified.  In  general,  the 
specification  of  available  potential  energy  and  available  internal  energy,  if 
present,  can  be  difficult,  particularly  for  those  models  that  include 
thermodynamic  effects  or  stratification.  For  further  information  see  Lorenz 
(1967)  and  Pedlosky  (1987)  among  others. 

Specification  of  the  energies  also  depends  on  the  volume  of  fluid  being 
analyzed.  For  example,  is  the  energy  analysis  to  be  made  for  a  unit  volume  of 
fluid,  a  column  of  fluid,  or  for  the  entire  region  being  studied?  For 
multilayer  hydrodynamic  models  it  seems  logical  to  integrate  kinetic  and 
available  potential  energy  densities  over  the  layer  depth  to  obtain  the 
kinetic  and  available  potential  energies  per  unit  area  for  each  layer.  This 
procedure  is  carried  out  in  sections  B.l  and  B.2  for  the  finite  depth  model 
and  is  called  the  layer-layer  formulation.  However,  as  shown  in  section  B.3, 
it  is  also  possible  to  associate  the  available  energy  with  a  given  interface; 
this  procedure  is  called  the  layer-interface  formulation.  Similarly,  for  the 
reduced-gravity  model,  available  potential  energy  per  unit  area  for  a  layer  or 
for  an  interface  may  be  used  as  shown  in  sections  B.5  and  B.6. 

The  second  problem  in  energy  analysis,  the  grouping  of  energy  transfer 
terms,  is  illustrated  for  the  finite  depth  model  by  comparing  the  kinetic 
energy  relations  for  the  layer-layer  formulation,  (2.27),  and  for  the  layer- 
interface  formulation,  (3.26).  Both  formulations  use  the  kinetic  energy  per 
unit  area  for  a  layer,  (1.1),  and  start  with  the  same  equation  for  the  local 
time  derivative  of  kinetic  energy  per  unit,  (1.11).  Thus,  both  the  kinetic 
energy  per  unit  area  in  a  layer  and  the  change  w.r.t.  time  of  this  quantity 
have  the  same  magnitude  in  the  two  formulations.  However,  the  layer-layer 
formulation  has  only  one  exchange  term  between  kinetic  and  available  potential 
energy  but  the  layer-interface  formulation  has  two. 

The  role  of  energetics  in  understanding  the  physical  mechanisms  present 
in  various  models  is  demonstrated  clearly  in  the  comparison  of  the  finite 
depth  and  reduced-gravity  models.  Consider  the  layer-layer  formulation  of  the 
energetics.  A  comparison  of  the  energy  block  diagrams  (Figs.  2  and  5)  shows 
that  the  same  type  of  energy  transfers  occur  in  both  models  but  a  comparison 
of  the  energy  relations,  (2.27)  -  (2.28)  and  (5.44)  -  (5.46),  shows  that  the 
magnitudes  of  these  transfer  terms  are  different.  The  situation  is  clarified 
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further  vhen  the  layer-interface  formulation  is  examined.  A  comparison  of  the 
energy  block  diagrams  (Figs.  3  and  6)  shows  that  the  energy  transfers  in  this 
formulation  are  different  for  the  two  models.  Specifically,  in  the  reduced- 
gravity  model  interface  available  potential  energy  is  transferred  between 
consecutive  interfaces,  but  no  such  transfer  occurs  in  the  finite  depth  model. 
This  direct  transfer  of  interface  available  potential  energy  between 
consecutive  interfaces  is  a  consequence  of  the  fact  that  the  pressure  in  the 
bottom  layer  of  a  reduced-gravity  model  must  be  independent  of  horizontal 
position  and  time.  This  independence,  in  turn,  dictates  that  the  interface 
deviation  at  the  top  of  the  bottom  layer,  Vu,  is  determined  completely  by  the 
density  structure  and  the  rest  of  the  interface  deviations. 
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Part  C,  Section  8 


The  energy  relations  in  the  layer-layer  formulation  are  given  in  this 
section  for  the  special  case  of  the  one-layer,  finite  depth  model.  A  vertical 
cross-section  of  the  one-layer  model  is  shovn  in  Figure  7. 

For  the  one-layer  model,  N=1  and  i=l.  From  (1.1)  and  (1.22),  the  kinetic  and 
available  potential  energies  are 

-  i  hjp,  (u/  .  v/)  (8.1) 

and 


Ai  -  2 


(8.2) 


Because  there  is  only  one  layer,  the  total  kinetic  energy  per  unit  area  and 
the  total  available  potential  energy  per  unit  area  are  represented  by  (8.1) 
and  (8.2),  respectively. 

From  (2.27)  and  (2.28),  for  i=l  and  N=l,  the  energy  relations  for  a  one  layer, 
finite  depth  model  are 


3Kj 
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3t  ■ 

Pi  ’ 

(2) 

where 

,  from  (2.5) 

and  Figure  7 

Pi  = 

Pig'll  =  P^g 

K  -  “i  ■] 

Term  1  is  the  energy  flux  divergence  and  term  2  is  the  exchange  between 
kinetic  energy  and  available  potential  energy  within  the  layer.  These 
relations  are  shown  in  the  block  diagram  in  Figure  8. 
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Adding  (8.3)  and  (8.4)  gives  the  local  time  derivative  of  the  total 
energy  per  unit  area  in  layer  1  as 

It  K  "  ''i]  =  -  ih  *  "ih]  • 

For  the  one  layer  model,  (8.6)  also  represents  the  local  time  derivative  of 
the  total  energy  per  unit  area  of  the  system.  As  expected  from  (2.32),  the 
only  term  on  the  RHS  is  the  energy  flux  divergence. 
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Part  C,  Section  9 


In  this  section  the  energy  relations  are  given  for  the  special  case  of  a 
two-layer,  finite  depth  model  in  the  layer-layer  formulation.  A  vertical 
cross-section  of  the  two-layer  model  is  shown  in  Figure  9. 


From  (1.1)  and  (1.22),  the  kinetic  and  available  potential  energies  per  unit 
area  in  a  layer  for  N=2  and  i=l,2  are 


^1  =  2  ^1^  K  "  ('^2  ~  ^*^2^l)]  ’ 
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(9.2) 

(9.3) 

(9.4) 


The  total  kinetic  energy  per  unit  area  for  the  two-layer,  finite  depth  model 
is  found  by  adding  (9.1)  and  (9.3)  to  get 
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(9.5) 


Similarly  the  total  available  potential  energy  per  unit  area  for  the  two- 
layer,  finite  depth  model  is  found  by  adding  (9.2)  and  (9.4)  to  find 
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From  (2.27)  and  (2.28),  for  i=l,2  and  N=2,  the  energy  relations  for  a  two- 
layer,  finite  depth  model  are 


fJjKi  .  hjPj) 

(1) 


^1  at  ’ 

(2) 


(9.7) 


ah^ 

ah^ 

at 

=  Pi  ar " 

(2) 

(3) 

(9.8) 


50 


(9.9) 
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where  the  pressures,  and  P2?  are  given  from  (2.5)  for  i=l,2  as 
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The  energy  transfers  in  (9.7)  -  (9.10)  are  shown  in  the  block  diagram  in 
Figure  10. 

Summing  (9.7)  and  (9.9)  gives  the  local  time  derivative  of  total  kinetic 
energy  per  unit  area  for  the  system  for  the  two-layer,  finite  depth  model  as 

It  K  ^2)  "  "  ^"^iK  ^1^1)  “  ^1  at  ^”^2(^2  ^2^2)  "  ^2  aF  ■ 
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Summing  (9.8)  and  (9.10)  gives  the  change  in  total  available  potential  energy 
per  area  for  the  two-layer  model  as 

/  V  9h.  ah„ 
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Note  that  (9.14)  has  no  term  3,  which  is  the  transfer  of  available  potential 
energy  between  the  two  layers.  Alternatively,  the  change  in  total  energy  per 
unit  area  in  a  layer  for  layers  1  and  2  can  be  found  by  adding  (9.7)  to  (9.8) 
and  (9.9)  to  (9.10),  respectively,  to  get 
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and 
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Here  the  exchanges  between  kinetic  and  available  potential  energy  with  a 
layer,  terms  (2)  and  (5),  are  missing.  The  change  with  time  in  total  energy 
per  unit  area  for  the  two-layer,  finite  depth  model  is  found  by  summing  (9.7) 
-  (9.10)  to  get 
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As  expected  from  (2.32), 
flux  divergences  for  the 


the  only  terms  on  the  RHS  of  (9.17)  are 
two  layers. 


the  energy 
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Part  C,  Section  10 


In  this  section  the  energy  relations  are  given  for  the  special  case  of  a 
three-layer  model  in  the  layer-layer  formulation.  A  vertical  cross-section  of 
the  three-layer  model  is  shown  in  Figure  11. 

From  (1.1)  and  (1.22),  the  kinetic  and  available  potential  energies  per 
area  in  a  layer  for  N=3  and  i=l,2,3,  are 
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(10.1) 
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The  total  kinetic  energy  per  unit  area  for  the  three-layer,  finite  depth  model 
is  found  by  adding  (10.1),  (10.3)  and  (10.5)  to  get 
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Similarly,  the  total  available  potential  energy  per  unit  area  for  the  three- 
layer,  finite  depth  model  is  found  by  adding  (10.2),  (10.4)  and  (10.6)  to  get 
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From  (2.27)  and  (2.28),  for  i=l,2,3  and  N=3,  the  energy  relations  for  a 
three-layer,  finite  depth  model  are 
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where  the  pressures,  p^,  p^,  are  given  from  (2.5)  for  i=l,2,3  as 
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These  energy  relations,  (10.9)  -  (10.14),  are  shown  in  the  block  diagram  in 
Figure  12. 


Summing  (10.9),  (10.11)  and  (10.13)  gives  the  change  in  the  total  kinetic 
energy  per  unit  area  of  the  system  for  the  three-layer,  finite  depth  model  as 
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(10.18) 
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Similarly,  summing  (10.10),  (10.12)  and  (10.14)  gives  the  change  in  the  total 
available  potential  energy  per  unit  area  for  the  three-layer  model  as 


3t  K  ^  ^  ^3) 
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(10.19) 


Alternatively,  the  change  in  total  energy  per  unit  area  in  a  layer  for  layers 
1-3  can  be  found  by  adding  (10.9)  to  (10.10),  (10.11)  to  (10.12),  and  (10.13) 
to  (10.14),  respectively,  to  get 


-  1 
3t  1 

h  * 

-  -  1 

S  *  ’’iPi)  " 

(piEhi] 

r3h2 

Ist  ^ 

9^3' 

yr. 

(10.20) 

(1) 

(3) 

^  1 
3t  1 

K  *  ''2] 

1  -  -  1 

[k2  .  h2P2)  t 

(piEhj  . 

3h3 

3t 

(4) 

(6) 

f  ^  1  (“'’2 

bt  *  at  , 

(10.21) 

(3) 

3 

3t 

(■^3  *  *3' 

)  =  - 

(S  ■"  ^3^3) 

(Plghl  . 

3h, 

1  ar 

3 

(10.22) 

(7) 

(6) 

The  change  in  the  total  energy  per  unit  area  for  the  three-layer  model  is 
found  by  summing  (10.9)  -  (10.14)  to  get 
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9t 


(kj  .  Aj  *  *  Kj  .  A3)  .  -  (kj  .  hjPj] 


(1) 


-  ^•V2  [k2  +  h^  P2)  ^•V3  [k^  +  p^j  . 

(4)  (7) 


Again,  as  expected  from  (2.32),  the  only  terms  on  the  RHS  of  (10.23) 
energy  flux  divergences  for  the  three  layers. 


(10.23) 


are  the 
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Part  C,  Section  11 


In  this  section  the  energy  relations  are  given  for  the  special  case  of  a 
one  active  layer,  reduced-gravity  model  in  the  layer-layer  formulations.  A 
vertical  cross-section  of  the  one  active  layer,  reduced-gravity  model  is  shown 
in  Figure  13. 

From  (1.22)  for  i=l,2  and  N=2 


^1  ■  2  [*^1  ('^2  ^*^2^l]. 

and 

A2  =  2  P2^  (*l2  ■  ^'^2^1]]  ■  ^ 

From  (1.1)  for  i=l,  and  from  the  discussion  following  (5.3),  the  kinetic 
energy  per  unit  area  in  layers  1  and  2  can  be  written  as 

Kj  .  I  hjPj  (u/  +  ( 


and 


K2  =  0  .  (11-M 

The  total  kinetic  energy  per  unit  area  for  the  one  active  layer,  reduced- 
gravity  model  is  found  by  adding  (11.3)  and  (11.4)  to  get 

Ki  *  K3  .  i  ,  (11.5, 

Similarly  the  total  available  potential  energy  per  unit  area  for  the  one 
active  layer,  reduced-gravity  model  is  found  by  adding  (11.1)  and  (11.2)  to 
get 

Ai  +  A2  =  ^  Pig  +\  (^2  “  ^1]^  ('^2^  “  ^'^2^1]  ■  (11*6) 

From  (5.44)  -  (5.46),  for  i=l,2  and  N=2,  the  energy  relations  for  a  one 
active  layer,  reduced-gravity  model  are 
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where,  from  (5.35)  and  (2.5),  for  i=l 


Pf  = 


and,  from  (5.22)  for  N=2 


(11.11) 


P 


2 


(11.12) 


The  energy  transfers  in  (11.7)  -  (11.10)  are  shown  in  the  block,  diagram  in 
Figure  14. 


Adding  (11.7)  and  (11.9)  gives 
unit  area  for  the  one  active  layer, 


the  change  in  the  total  kinetic  energy  per 
reduced-gravity  model  as 


h  K  "  ''2) 


(kj  .  hjfj) 

(1) 


^1  at 
(2) 


(11.13) 


Similarly,  summing  (11.8)  and  (11.10)  gives  the  change  in  total  available 
potential  energy  per  unit  area  for  the  one  active  layer,  reduced-gravity  model 
as 
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Alternately,  the  change  in  total  energy  per  unit  area  in  layers  1  and  2  can  be 
found  by  adding  (11.7)  to  (11.8)  and  (11.9)  to  (11.10),  respectively,  to  get 
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Finally,  the  change  in  the  total  energy  per  unit  area  for  the  one  active 
layer,  reduced-gravity  model  is  found  by  summing  (11.7)  -  (11.10)  to  get 


It  K  +  +  A2]  -  -  ^-v^  +  h^pj] 

or,  from  (5.52)  for  N=2, 

It  K  +  +  K2  +  A^J  =  -  +  h^p^J  +  P2  JT  ' 


(11.17) 


(11.18) 
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Part  C,  Section  12 


In  this  section  the  energy  relations  are  given  for  the  special  case  of  a 
two  active  layer,  reduced-gravity  model  in  the  layer-layer  formulation.  A 
vertical  cross  section  of  the  two  active  layer,  reduced-gravity  model  is  shown 
in  Figure  15. 


From  (1.22)  for  i=l,2,3  and  N=3,  the  available  potential  energy  per  unit 
area  for  each  layer  is 


^1  "  2  hi  "^2  "  ^*^2^1]]  ’ 

^2^2  ^2^  [^2  "  ^'^2^1)  "  ^3  ■  ^*^3  K  ^2]].  ’ 
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(12.2) 

(12.3) 


From  (1.1)  for  i=l,2,  and  from  the  discussion  following  (5.3),  the  kinetic 
energy  per  unit  area  in  layers  1,2,  and  3  can  be  written  as 


^2  =  2  "2^2  (^2  ^"2  ]  ’ 


(12.4) 
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(12.6) 


The  total  kinetic  energy  per  unit  area  for  the  two  active  layer,  reduced- 
gravity  model  is  found  by  summing  (12.4)  -  (12.6)  to  get 
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(12.7) 


Similarly,  the  total  available  potential  energy  per  unit  area  for  the  two 
active  layer,  reduced-gravity  model  is  found  by  summing  (12.1)  -  (12.3)  to  get 


(12.8) 
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From  (5.44)  -  (5.46)  for  i=l,2,3  and  N=3,  the  energy  relations  for  a  two 
active  layer,  reduced-gravity  model  are 
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where,  from  (5.35)  and  (2.5),  for  i=l,2 
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and,  from  (5.20),  N=3 
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The  energy  relations,  (12.9)  -  (12.14),  are  shown  in  the  block  diagram  in 
Figure  16. 

Summing  (12.9),  (12.11),  and  (12.13)  gives  the  change  in  the  total 
kinetic  energy  per  unit  area  for  the  two  active  layer,  reduced-gravity  model 
as 
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Similarly,  summing  (12.10),  (12.12)  and  (12.14)  gives  the  change  in  the  total 
available  potential  energy  per  unit  area  for  the  two  active  layer,  reduced- 
gravity  model  as 
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Alternatively,  the  change  in  total  energy  per  unit  area  in  layers  1,2,  and  3 
can  be  found  by  adding  (12.9)  to  (12.10),  (12.11)  to  (12.12),  and  (12.13)  to 
(12.14),  respectively,  to  get 


(12.20) 
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Finally,  the  change  in  the  total  energy  per  unit  area  for  the  two  active 
layer,  reduced-gravity  model  is  found  by  summing  (12.9)  -  (12.14)  to  get 
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or,  from  (5.52)  for  N=3, 
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Part  D,  Section  13 


In  this  section  the  energy  specifications  and  relations  are  given  for  the 
special  case  of  a  one-layer,  finite  depth,  model  in  the  layer-interface 
formulation.  The  reader  should  be  familiar  with  section  B.3  where  the  layer- 
interface  formulation  is  discussed  in  detail  for  the  finite  depth  model.  A 
vertical  cross-section  of  the  one-layer  model  is  shown  in  Figure  7.  From 
(1.1),  for  i=l  and  N=l,  the  kinetic  energy  per  unit  area  for  the  layer  is 

Kj  -  I  hjpj  .  v/)  .  (13.1) 

For  the  interface  formulation,  the  available  potential  energy  per  unit  area  is 
found  in  (3.4).  For  i=l  and  N=l,  this  relation  is 
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From  (3.26)  and  (3.27),  for  i=l  and  N=l,  the  energy  relations  for  the  one 
layer  model  in  layer-interface  formulation  are 
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where,  from  (2.5)  and  Figure  7, 

Pi  =  Pighi  =  Pig  (hi  -  hJ  . 

The  energy  transfers  in  (13.3)-(13.4)  are  shown  in  the  block,  diagram  in  Figure 
17.  Adding  (13.3)  and  (13.4)  gives 
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as  expected  from  (3.33)  for  N=l. 

Eliminate  from  (13.3)  and  (13.4)  by  using  (13.5), 
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Note  that  the  energy  relations  for  the  one  layer  model  in  the  layer-layer 
formulation,  (8.3)  -  (8.4),  are  identical  to  those  in  the  layer-interface 
formulation,  (13.7)  -  (13.8).  Furthermore,  the  total  kinetic  energy,  (8.1) 
and  (13.1),  and  the  total  available  potential  energy,  (8.2)  and  (13.2),  are 
identical  for  the  two  formulations  as  expected. 
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Part  D,  Section  14 


In  this  section  the  energy  relations  are  given  for  the  special  case  of  a 
two-layer,  finite  depth  model  in  the  layer-interface  formulation.  A  vertical 
cross  section  of  the  two  layer  model  is  shown  in  Figure  9.  From  (1.1),  for 
i=l,2  and  N=2,  the  kinetic  energy  per  unit  area  in  layers  1  and  2  can  be 
written  as 

Kj  -  I  hjp,  (U.l) 

and 

K2  =  I  h2P2  (u/  +  V2^)  (1^-2) 

From  (3.4),  the  available  potential  energy  per  unit  area  in  the  interface 
formulation  for  layers  1  and  2  is  expressed  as 


-  2  >^1 


and 


(14.3) 


h  =  H'’2  -  Pi)®  (^2^  -  ^^2^1)  ’ 


(14.4) 


The  total  kinetic  energy  per  unit  area  for  this  case  is  found  by  adding  (14.1) 
and  (14.2)  to  get 


K  ^2)  "  2  ^1^1  K  ^^1)^2  ^2^2  (“2  ^  ^2  ] 


(14.5) 


Similarly  the  total  available  potential  energy  per  unit  area  for  this  case  is 
found  by  adding  (14.3)  and  (14.4)  to  find 


-  -  1  2  1  f 

A)  +  A2  =  2  Pi^  >^1  +  2  r2 


(14.6) 


Note  that  the  total  kinetic  energy  for  the  two-layer  case  is  the  same  in  the 
layer-layer  formulation,  (9.5),  and  in  the  layer-interface  formulation.  The 
total  available  potential  energy  is  the  same  in  both  formulations,  (9.6)  vs 
(14.6). 

From  (3.26)  and  (3.27),  for  i=l,2  and  N=2,  the  energy  relations  for  a 
two-layer  model  in  the  layer-interface  formulation  are 
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where  the  pressures,  and  P2>  are  given  from  (2.5)  with  i=l,2  as 

(14.11) 

h  =  -  (P2  - 

The  energy  transfers  in  (14.7)  -  (14.10),  are  shown  in  Figure  18. 


Summing  (14.7)  and  (14.9)  gives  the  change  in  total  available  potential 
energy  per  unit  area  for  this  case  that  is 
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From  Figure  9, 
depth,  h^,  for 


the  relations  between  interface  deviation, 
i=l  and  2  are  expressed  as 


0^,  and  the  layer 


hj^  =  -  ri2 


(14.14) 


and 
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—  H2  +  ^2  ■ 


(14.15) 


Substituting  (14.14)  -  (14.15)  into  (14.13)  yields 


(14.16) 


A  comparison  of  (14.16)  and  (9.14)  reveals  that  the  change  in  total  available 
potential  energy  per  unit  area  for  the  two-layer  case  is  the  same  for  both  the 
layer-layer  and  the  layer-interface  formulations. 

The  energetics  for  a  two-layer,  finite  depth  model  have  been  derived  by 
Hurlburt  and  Thompson  (1982).  They  specify  the  kinetic  energy  per  unit  area 
for  each  layer,  but  use  the  total  available  potential  energy  per  unit  area  for 
the  system.  Thus,  the  question  of  which  formulation  to  use  for  the  available 
potential  energy  does  not  arise.  The  relations  derived  here  are  equivalent  to 
those  in  their  paper  with  three  differences:  (1)  They  use  p  instead  of 
and  when  these  values  appear  alone,  (2)  the  relations  derived  here  have  no 
source  or  sink  terms  included,  and  (3)  the  relations  derived  here  ((14.8), 
(14.10),  and  (14.13))  contain  extra  terms  that  vanish  if  averaged  over  the 
entire  model  region. 

Holland  and  Lin  (1975)  also  have  examined  the  energetics  for  a  two-layer, 
finite  depth  model.  They  specify  the  kinetic  energy  per  unit  area  for  each 
layer  and  the  total  available  potential  energy  per  unit  area  for  the  system. 
Then  they  average  these  quantities  over  the  entire  basin.  However,  they  use  a 
rigid  top  (hi  =  0)  with  their  model  so  the  area-averaged,  total  available 
potential  energy  is  simply  the  area-averaged,  available  potential  energy 
associated  with  the  interface  between  the  two  layers. 

Summing  (14.8)  and  (14.10)  yields  the  change  in  total  kinetic  energy  per 
unit  area  for  this  case  as 
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Substituting  (14.14)  -  (14.15)  into  (14.17)  gives 
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(14.18) 
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which  is  the  same  as  (9.13),  the  change  in  total  kinetic  energy  per  unit  area 
in  the  layer-layer  formulation  for  the  two-layer  case. 

The  change  in  total  energy  per  unit  area  for  the  system  is  found  by  summing 
(14.7)  -  (14.10)  to  get 
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(2)  (5) 


(14.19) 


as  expected  from  (3.33).  Note  that  (14.19)  is  identical  to  (9.17),  which 
shows  that  the  change  in  total  energy  per  unit  area  for  the  two-layer  case  of 
the  finite  depth  model  does  not  depend  on  the  type  of  formulation  for 
available  potential  energy. 
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Part  D,  Section  15 


In  this  section  the  energy  specifications  and  relations  are  given  for  the 
special  case  of  a  three-layer,  finite  depth  model  in  the  layer-interface 
formulation.  A  vertical  cross  section  of  the  three-layer  model  is  shown  in 
Figure  11.  From  (1.1),  for  i=l,2,3  and  N=3,  the  kinetic  energy  per  unit  area 
in  layers  1,  2,  and  3  can  be  written  as 
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1  h^Pl  , 

2  ^3^3  (“3  ^  ^"3  )  • 


(15.1) 

(15.2) 

(15.3) 


The  total  kinetic  energy  per  unit  area  for  this  case  is  found  by  summing 
(15.1)  -  (15.3)  to  get 
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From  (3.4),  for  i=l,2,  and  3  and  N=3,  the  available  potential  energy  per 
unit  area  for  layers  1,  2  and  3  can  be  expressed  as 
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The  total  available  potential  energy  per  unit  area  for  this  case  is  found  by 
summing  (15.5)  -  (15.7)  to  get 
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Note  that  the  total  kinetic  energy  is  the  same  for  the  three-layer  model  in 
the  layer-layer  formulation,  (10.7),  and  the  layer-interface  formulation, 
(15.4).  Similarly,  the  total  available  potential  energy  per  unit  area  is  the 
same  for  both  formulations,  (10.8)  and  (15.8). 

From  (3.26)  and  (3.27),  for  i=l,  2  and  3  and  N=3,  the  energy  relations 
for  a  three-layer,  finite  depth  model  in  the  layer-interface  formulation  are 
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where 

the  pressure,  p^ , 

is  given  in 

(2.5)  for 

II 

(15.9) 


(15.10) 


(15.11) 


(15.12) 


(15.13) 
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(15.15) 

(15.16) 
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and 


P3  =  P3g>l3  -  (P3  -  P 


ijg^i  -  (p3  P2)sh2  • 


(15.17) 


These  energy  relations,  (15.9)  -  (15.14),  are  shown  in  the  block  diagram  in 
Figure  19. 

Summing  (15.9),  (15.11)  and  (15.13)  gives  the  change  in  the  total 
available  potential  energy  per  unit  area  for  this  case  that  is 
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From  Figure  11,  the  relations  between  interface  deviation,  and  the  layer 


depth,  h^,  for  i=l,2,  and  3  are  expressed  as 


h2  ~  ^2  ^2  ^3  ^ 


^3  “  *^3  ^3  ■ 


(15.19) 

(15.20) 

(15.21) 


Substituting  (15.19)  -  (15.21)  into  (15.18)  yields 

3h,  ah„  8h, 


+  A2  +  A^ 


)  .  '’I  .  2  .  3 

J  "  Pi  at  P2  at  P3  at  ' 


(15.22) 


A  comparison  of  (15.22)  and  (10.19)  reveals  that  the  change  in  total  available 
potential  energy  per  unit  area  for  the  three-layer  case  of  the  finite  depth 
model  is  the  same  for  the  layer-layer  and  layer-interface  formulations. 

The  change  w.r.t.  time  in  the  total  kinetic  energy  per  unit  area  is  found 
by  adding  (15.10),  (15.12),  and  (15.14).  The  result  is 
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(15.23) 
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Substituting  (15.19)  -  (15.21)  into  (15.23)  yields 

3h, 


It  K  *  >^2  *  '^3)  ■  -  K  *  '■iPi)  ■  h  JT 

-  |k2  +  h2P2j  -  P2 

-  ^-^3  [h  ^  *^3^]  -  hJT  ’ 


(15.24) 


which  is  the  same  as  (10.18),  the  change  in  total  kinetic  energy  per  unit  area 
in  the  layer-layer  formulation  for  the  three-layer,  finite  depth  model. 

The  change  in  total  energy  per  unit  area  for  the  system  is  found  by 
summing  (15.9)  -  (15.14)  to  get 

(Aj  .  Kj  .  Aj  .  *  *3  *  *"3]  ■  -  *-“1  h  *  '’if'l) 

(2) 

-  ^•V2  (k2  +  h2P2)  -  ^-^3  [K3  +  h^p^]  ,  (15.25) 

(5)  (8) 

as  expected  from  (3.33)  for  N=3.  Note  that  (15.25)  is  identical  to  (10.23), 
which  shows  that  the  change  in  total  energy  per  unit  area  for  the  three-layer, 
finite  depth  model  does  not  depend  on  the  type  of  formulation. 
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Part  D,  Section  16 


In  this  section,  the  energy  specifications  and  relations  are  given  for 
the  special  case  of  a  one  active  layer,  reduced-gravity  model  in  the  layer- 
interface  formulation.  The  reader  should  be  familiar  with  section  B.6  where 
the  layer-interface  formulation  for  the  reduced-gravity  model  is  discussed  in 
detail.  A  vertical  cross-section  of  the  one  active  layer,  reduced-gravity 
model  is  shown  in  Figure  13.  From  (1.1),  for  i=l,  and  from  the  discussion 
following  (5.3),  the  kinetic  energy  per  unit  area  in  layers  1  and  2  can  be 
written  as 

^1  =  5  [u^^  +  Vj^j  (16.1) 


and 


K 


2 


=  0  . 


(16.2) 


The  total  kinetic  energy  per  unit  area  for  the  layers  in  the  reduced-gravity 
model  is  found  by  adding  (16.1)  and  (16.2)  to  get- 


V  V  Ik  r  2  2) 

^1  ^2  ~  2  ^1^1  t  1  '^l  j 


From  (3.4)  for  i=l,2  and  N=2. 


Af  -  2  hj 


^2  =  2  (^2  Pl]^  ('^2  ^'^2^l) 


(16.3) 


(16.4) 

(16.5) 


The  total  available  potential  energy  per  unit  area  for  this  case  is  found  by 
adding  (16.4)  and  (16.5)  to  obtain 


^1  ^2  ~  2  ^1®  ^1  ^  2  1^2  ^1^^  (^2  ~  ^^2^lj 


(16.6) 


Note  that  the  total  kinetic  energy  per  unit  area  is  the  same  for  the  one 
active  layer,  reduced-gravity  model  in  the  layer-layer  formulation,  (11.5), 
and  the  layer-interface  formulation,  (16.3).  The  total  available  potential 
energy  per  unit  area  for  this  case  is  the  same  for  both  formulations,  (11.6) 
and  (16.6). 

From  (6.21)  -  (6.23),  for  i=l,2  and  N=2,  the  energy  relations  for  a  one 
active  layer,  reduced-gravity  model  in  the  layer-interface  formulation  are 
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where,  from  (6.11)  and  (2.5),  for  i=l, 


Pi  - 

P2  =  Pi  - 

and,  from  (6.4)  and  N=2, 

P2  =  (^1  “  ‘’2)®^! 


(16.7) 


(16.8) 


(16.9) 


(16.10) 


(16.11) 


(16.12) 


These  energy  relations,  (16.7)  -  (16.10),  are  shown  in  the  block,  diagram  in 
Figure  20. 


Summing  (16.8)  and  (16.10)  gives  the  change  in  the  total  kinetic  energy 
per  unit  area  for  the  one  active  layer,  reduced-gravity  model  in  the  layer- 
interface  formulation  as 


It  h  *  ‘'2] 


f?!  (Kj  .  hjfj) 

(3) 


3n^  ^  302 
Pi  ^  Pi  3F 
(1)  (4) 


(16.13) 


From  Figure  13,  the  relation  between  interface  deviation,  n^j  and  depth,  h^, 
for  i=l,2  is 
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-  ri2  • 


(16.14) 


Using  (16.14)  in  (16.13)  yields 

It  Vl  ^2]  "  “  ^*^^1  i^l  ^1^1]  "  ^1  It” 


(16.15) 


which  is  the  same  as  (11.13),  the  change  in  total  kinetic  energy  per  unit  area 
in  the  one  active  layer,  reduced-gravity  model  in  the  layer-layer  formulation. 

Summing  (16.7)  and  (16.9)  gives  the  change  in  the  total  available 
potential  energy  per  unit  area  for  this  case,  which  is 
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Substituting  (16.14)  into  (16.16)  yields 
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A  comparison  of  (16.17)  and  (11.14)  shows  that  the  change  in  total  available 
potential  energy  per  unit  area  for  the  one  active  layer,  reduced-gravity  model 
does  not  differ  between  the  layer-layer  and  the  layer-interface  formulations. 


The  change  in  the  total  energy  per  unit  area  for  this  case  is  found  by 
summing  (16.7)  -  (16.10)  to  find 


It  (^1  +  +  ^2  +  K2]  -  -  ^-v^  (k^  +  h^pj 
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(16.18) 


or,  from  (6.28) 
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(16.19) 


Note  that  (16.19)  and  (11.18)  are  identical  so  that  the  change  in  total  energy 
per  unit  area  does  not  depend  on  the  type  of  formulation:  layer-layer  or 
layer-interface . 
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Part  C,  Section  17 


In  this  section,  the  energy  specifications  and  relations  are  given  for 
the  special  case  of  a  two  active  layer,  reduced-gravity  model  in  layer- 
interface  formulation.  A  vertical  cross-section  of  the  two  active  layer, 
reduced-gravity  model  is  shown  in  Figure  15.  From  (1.1)  for  i=l,2  and  from 
the  discussion  following  (5.3)  for  N=3,  the  kinetic  energy  per  unit  area  in 
layers  1,2  and  3  can  be  written  as 


Kl  .  1  hjPj  , 

*^2  ■  2  ''2  ”2  (“2  *  ''2  )  ’ 


(17.1) 

(17.2) 


(17.3) 


The  total  kinetic  energy  per  unit  area  for  all  layers  in  this  reduced-gravity 
model  is  found  by  summing  (17.1)  -  (17.3)  to  get 
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(17.4) 


From  (3.4)  for  i=l,2,3  and  N=3,  the  available  potential  energy  in  the 
interface  formulation  for  layers  1-3  can  be  written  as 


Ai  -  2  P^g 


(17.5) 


A2  -  2  (p2  -  Pi)e  ['^2  ^'^2^1)  ’ 

^3  "  2  (‘^3  “  ^2]^  [*^3  "  ^^^3  K  ^2].  ’ 


The  total  available  potential  energy  for  this  case  is  found  by  summing 
(17.5)  -  (17.7)  to  obtain 

^1  ^2  ^  ^3  ”  2  ^1^  ^1  ^2  (^2  ~  ^ij®  (^2  ”  ^^2^1) 

.  I  (pj  -  Pjjg  [^3=^  -  2113  (Hj  .  Hj]]  .  (17. 

Note  that  the  total  kinetic  energy  per  unit  area  is  the  same  for  the  two 
active  layer,  reduced-gravity  model  in  the  layer-layer  formulation,  (12.7), 
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and  for  the  layer-interface  formulation,  (17.4).  The  total  available 
potential  energy  per  unit  area  for  this  case  is  the  same  for  both 
formulations,  (12.8)  and  (17.8). 

From  (6.21)  -  (6.23)  for  N=3,  the  energy  relations  for  a  two  active 
layer,  reduced-gravity  model  in  the  layer-interface  formulation  are 
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where,  from  (6.11)  and  (2.5),  for  i=l  and  2, 
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(17.16) 


and,  from  (6.4),  for  N=3, 

P3  ■  (Pl  “  '3)®"!  *  (P2  -  P3)®“2  • 


(17.17) 


These  energy  transfers  in  (17 . 9)~(17 • 14)  are  shown  in  the  block  diagram  in 
Figure  21, 

Summing  (17-10),  (17.12),  and  (17.14)  gives  the  change  in  the  total 
kinetic  energy  per  unit  area  for  the  two  active  layer,  reduced-gravity  model 
in  the  layer-interface  formulation  as 
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From  Figure  15,  the  relation  between  interface  deviation,  »  and  layer  depth, 

h.,  can  be  written  for  layers  1  and  2  as 
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and 
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Using  (17.19)  -  (17.20)  in  (17.18)  yields 


(17.20) 
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3t 
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(17.21) 


which  is  the  same  as  (12.18),  the  change  in  total  kinetic  energy  per  unit  area 
for  the  two  active  layer,  reduced-gravity  model  in  the  layer-layer 
formulation. 

Summing  (17.9),  (17.11)  and  (17.13)  gives  the  change  in  the  total  available 
potential  energy  per  unit  area  for  this  case,  which  is 
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Substituting  (17.19)  -  (17.20)  into  (17.22)  yields 
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(17.23) 


A  comparison  of  (17.23)  and  (12.19)  shows  that  the  change  in  total  available 
potential  energy  for  the  two  active  layer,  reduced-gravity  model  does  not 
differ  between  the  layer-layer  and  the  layer-interface  formulations. 

The  change  in  the  total  energy  per  unit  area  for  this  case  is  found  by 
summing  (17.9)  -  (17.14)  to  obtain 
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Note  that  (17.24)  and  (12.23)  are  identical  as  are  (17.25)  and  (12.24)  so  that 
the  change  in  total  energy  per  unit  area  does  not  depend  on  the  type  of 
formulation:  layer-layer  or  layer-interface. 
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PART  E:  Summary  and  Conclusions 


Kinetic  energy  and  potential  energy  are  included  in  the  energy  analysis 
of  hydrodynamic  models.  Internal  energy  is  omitted  because  these  models  do 
not  use  temperature  as  a  variable.  The  energy  analysis  includes  the 
specification  of  energies  present  and  the  change  in  magnitude  of  these 
energies  with  v.r.t.  time.  Consider  a  fluid  system  that  is  surrounded  by 
solid  vails  and  bottom.  If  there  are  no  external  sources  and  sinks  of  energy, 
the  total  of  all  energies  present  in  the  system  remains  constant  with  time. 
However,  energy  may  be  redistributed  horizontally  and  vertically  within  the 
system  and  may  be  transferred  from  one  type  to  another.  The  ways  in  which 
energy  is  redistributed  and  transferred  give  insight  into  the  behavior  of  the 
dynamical  model. 

Usually  less  than  one  percent  of  the  total  potential  energy  can  be 
transferred  into  other  energy  types  or  redistributed  in  space,  so  this 
available  potential  energy,  rather  than  the  total  potential  energy,  is 
considered  in  the  energy  analysis  of  a  fluid  system.  This  available  potential 
energy  is  defined  as  the  difference  between  the  total  potential  energy 
actually  present  and  some  minimum  potential  energy  state.  Choosing  the 
minimum  potential  energy  state  can  be  difficult  for  some  dynamical  models. 

For  the  finite  depth  and  reduced-gravity  models  studied  here,  the  minimum 
energy  state  is  defined  as  the  state  with  no  motion  and  no  interface 
deviations. 

Both  kinetic  and  available  potential  energy  can  be  specified  per  unit 
volume,  but  for  layered  models  it  seems  reasonable  to  integrate  over  a  layer 
to  find  both  kinetic  and  available  potential  energy  per  unit  area  in  a  layer. 
This  layer  specification  of  energies  is  used  first  in  section  B.l,  and  is 
called  the  layer-layer  formulation.  However,  as  shown  in  section  B.3,  it  is 
also  possible  to  associate  available  potential  energy  with  a  given  interface. 
The  specification  of  the  kinetic  energy  per  unit  area  in  a  layer  and  the 
interface  available  potential  energy  is  called  the  layer-interface 
formulation.  The  vertical  sorting  of  available  potential  energy  is  different 
for  the  two  formulations,  but  the  total  available  potential  energy  per  unit 
area,  found  by  summing  over  the  vertical  extent  of  the  system,  is  identical 
for  both  formulations  as  it  should  be. 

Once  the  energies  are  specified,  the  change  with  w.r.t.  time  of  these 
energies  can  be  found  readily  from  the  basic  equations  as  shown  in  section  B.l 
for  the  finite  depth  model.  Then  the  terms  in  the  energy  relations  must  be 
grouped  carefully  to  give  the  energy  transfers  possible.  As  shown  in  section 
B.4,  it  is  possible  to  obtain  a  set  of  energy  transfers  that  are 
mathematically  correct,  but  physically  meaningless. 

For  the  finite  depth  model,  the  layer-layer  and  the  layer-interface 
formulations  both  lead  to  physically  reasonable,,  but  different,  energy 
transfers  as  shown  in  Figures  2  and  3.  For  the  layer-layer  formulation, 
transfer  between  kinetic  and  available  potential  energies  occurs  within  a 
layer,  but  only  available  potential  energy  is  exchanged  between  adjacent 
layers.  For  the  layer-interface  formulation,  kinetic  energy  in  a  layer  is 
exchanged  with  interface  available  potential  energy  for  the  interfaces 
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directly  above  and  below  the  layer.  The  interface  available  potential  energy 
is  exchanged  with  the  kinetic  energy  in  the  layers  directly  above  and  below 
the  interface.  The  magnitudes  of  both  the  kinetic  energy  per  unit  area  in  a 
layer  and  its  change  w.r.t.  time  are  identical  for  the  two  formulations.  Both 
formulations  offer  useful  insights  into  the  energetics  of  the  finite  depth 
model,  and  there  appears  to  be  no  reason  to  prefer  one  over  the  other. 

The  basic  equations  for  the  finite  depth  model,  together  with  the 
assumptions  and  approximations  are  found  in  the  Appendix.  The  reduced-gravity 
model,  discussed  in  detail  in  section  B.5,  has  the  additional  assumption  that 
the  bottom  layer  is  infinitely  deep.  The  kinetic  energy  per  unit  area  for  an 
infinitely  deep  layer  is  finite  only  if  the  velocity  in  the  layer  is  zero. 

The  condition  of  zero  velocity,  together  with  conservation  of  mass,  means  that 
the  pressure  in  the  bottom  layer  does  not  vary  with  time  or  horizontal 
position.  Therefore,  the  interface  deviations  cannot  all  be  specified 
independently  or,  in  other  words,  the  interface  deviation  at  the  top  of  the 
bottom  layer  is  determined  by  the  density  structure  and  the  rest  of  the 
interface  deviations. 

In  the  layer-layer  formulation,  the  same  energy  transfers  are  present  for 
the  finite  depth  and  reduced-gravity  models  (Figs.  2  and  5),  but  the 
magnitudes  of  these  transfers,  particularly  for  the  exchange  of  available 
potential  energy  between  adjacent  layers,  are  different.  In  the  layer- 
interface  formulation,  the  energy  transfers  for  the  two  models  (Figs.  3  and  6) 
differ  in  type  as  well  as  magnitude.  For  the  reduced-gravity  model  only, 
interface  available  potential  energy  is  exchanged  between  the  interfaces  above 
and  below  a  layer.  At  first  it  may  appear  physically  unsatisfactory  for 
interface  available  potential  energy  to  be  exchanged  between  two  interfaces 
that  are  not  in  physical  contact.  However,  this  exchange  merely  reflects  the 
fact  that  in  the  reduced-gravity  model  the  interface  deviations  and, 
therefore,  the  interface  available  potential  energies,  are  coupled  and  are  not 
all  independent.  Thus,  both  formulations  are  reasonable  physically,  but  the 
greater  complexity  of  the  layer-interface  formulation  might  make  the  layer- 
layer  formulation  more  desirable  for  the  reduced-gravity  model. 

In  conclusion,  two  different,,  but  physically  consistent,  formulations  of 
the  energetics  for  the  finite  depth  and  the  reduced-gravity  models  have  been 
found. 
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Appendix:  Derivation  of  the  basic  equations 


Consider  the  equations  for  conservation  of  mass  and  momentum  in  a 
multilayer  hydrodynamic  ocean  model.  Assume  the  fluid  is  inviscid, 
hydrostatic,  and  rotating.  In  addition,  assume  that  the  density  in  each  layer 
is  held  constant.  The  resulting  equations  for  conservation  of  mass  and 
momentum  in  a  given  layer,  i,  are 
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(A. 2) 


(A. 3) 


(A. 4) 


Subscripts  indicate  the  layer;  other  symbols  are  defined  in  the  Glossary. 
Figure  1  is  a  vertical  cross  section  of  the  multilayer  model.  The  pressure  at 
depth  z  can  be  obtained  by  integrating  (A. 3)  vertically  from  depth  z  to  the 
top  of  the  ocean: 

\  ^ 

Pjgdz  (A. 5) 

2 


or 


p.(x,y,z,t)  =  p^(x,y,t)  +  Pig(n^  -  zj  -  Z  (p.  -  pjgh^  (A. 6) 


The  applied  atmospheric  pressure,  p  (x,y,t),  will  be  set  to  zero  and  ignored 
for  these  calculations.  Note  that  The  only  z  dependence  on  the  (RHS)  of  (A. 6) 
occurs  in  the  p^gz  term.  Thus,  (A. 6)  may  be  written  as 

Pj(x,y,z,t)  =  p^(x,y,t)  -  p.gz  ,  (A. 7) 

where 
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Pi(x,y,t)  =  p.grij  -  [p.  -  pj]ghj 


(A. 8) 


Operating  on  (A. 7)  with  where  ^  ^  +  j  yields 


^  p.(x,y,z,t)  =  ^  P,.(x,y,t)  . 


(A. 9) 


Thus,  there  is  no  z  dependence  on  the  RHS  of  (A.l)  and  (A. 2).  It  follows  that 
if  the  horizontal  velocity  is  independent  of  z  initially,  it  will  remain  that 
way.  Then  (A.l)  and  (A. 2)  become 
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Integrating  equation  (A. 4)  over  the  thickness  of  a  layer  yields 
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(A. 12) 


The  horizontal  velocity  does  not  vary  with  z  within  a  layer,  so  (A. 12)  becomes 
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Using  (A. 14)  and  (A. 15)  in  (A. 13)  yields 


86 


(A. 16) 


(A. 6),  (A. 10),  (A. 11)  and  (A. 17)  form  the  set  of  four  equations  and  four 
unknowns  for  each  layer  of  the  various  layered  ocean  models  examined  in  this 
study.  This  set  will  be  used  to  derive  the  energy  transformations  for  each 
model  and  is  listed  for  easy  reference  in  Table  1  on  page  4. 
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Figure  1.  A  vertical  cross  section  of  an  N-layer,  finite  depth  model. 
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Figure  2.  Energy  transfers  for  three  adjacent  layers  in  an  N-layer,  finite 
depth  model  for  the  layer-layer  formulation. 
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Figure  3.  Energy  transfers  for  three  adjacent  layers  and  interfaces  in 
an  N-layer,  finite  depth  model  for  the  layer-interface  formulation. 
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Figure  4.  Energy  transfers  for  three  adjacent  layers  and  interfaces  in  an 
N-layer,  finite  depth  model  for  a  physically  unrealistic  formulation. 
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Figure  5.  Energy  transfers  for  the  bottom  three  layers  in  an  N-layer, 
reduced-gravity  model  for  the  layer-layer  formulation. 
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Figure  6.  Energy  transfers  for  the  bottom  three  layers  and  interfaces  in  an 
N~layer,  reduced-gravity  model  for  the  layer-interface  formulation* 
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Figure  7.  A  vertical  cross-section  of  a  one-layer,  finite  depth  model. 


Figure  8.  Energy  transfers  in  a  one-layer,  finite  depth  model  for  the  layer- 
layer  formulation. 
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Figure  9.  A  vertical  cross-section  of  a  two-layer,  finite  depth  model. 
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Figure  10.  Energy  transfers  in  a  two-layer,  finite  depth  model  for  the  layer- 
layer  formulation. 
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Figure  11.  A  vertical  cross  section  of  a  three-layer,  finite  depth  model. 
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Figure  12.  Energy  transfers  in  a  three-layer,  finite  depth  model  for  the 
layer-layer  formulation. 
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Figure  13.  A'vertical  cross  section  of  a  one  active  layer,  reduced-gravity 
inod6l  • 
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Figure  14.  Energy  transfers  in  a  one  active  layer,  reduced-gravity  model  for 
the  layer-layer  formulation. 


Figure  15.  A  vertical  cross  section  of  a  two  active  layer,  reduced-gravity 
model. 


97 


Fi 


A. 

-  ©  1 

- Z - 

Ao 

K2 

0  ^ 

— I — 

6 


Figure  16.  Energy  transfers  in  a  two  active  layer,  reduced-gravity  model  for 
the  layer-layer  formulation. 


Figure  17.  Energy  transfers  in  a  one-layer,  finite  depth  model  for  the  layer- 
interface  formulation. 
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Figure  18.  Energy  transfers  in  a  two-layer,  finite  depth  model  for  the  layer- 
interface  formulation. 
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Figure  19.  Energy  transfers  in  a  three-layer,  finite  depth  model  for  the 
layer-interface  formulation. 
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Figure  20.  Energy  transfers  in  a  one  active  layer,  reduced-gravity  model  for 
the  layer-interface  formulation. 
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Figure  21.  Energy  transfers  in  a  two  active  layer,  reduced-gravity  model  for 
the  layer-interface  formulation. 
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